Chinav :人 
C hinaA IVE 


心理 科学 进展 2021, Vol. 29, No. 11, 1948-1969 © 2021 中 国 科学 院 心理 研究 所 
Advances in Psychological Science https://doi.org/10.3724/SP.J.1042.2021.01948 


* 9t i, a ik (Research Method) * 


密集 追踪 数据 分 析 : 模型 及 其 应 用 


MEH KHAS FRE BRR 
(中 山大 学 心理 学 系 , 广州 510006) 


摘 要 在 心理 学 、 教 育 学 和 临床 医学 等 领域 ， 越 来 越 多 的 研究 者 开始 关注 个 体内 部 的 行为 、 心 理 、 临 床 效 
果 等 随时 间 而 产生 的 动态 变化 ,重视 针对 个 体 的 差异 化 建 模 。 密 集 追 踪 是 一 种 在 短 时 间 内 对 个 体 进 行 多 个 时 
间 节 点 密集 追踪 测量 的 方法 ， 更 适合 用 于 研究 个 体内 部 心理 过 程 等 的 动态 变化 及 其 作用 机 制 。 近 年 来 ， 密 集 追 
踪 成 为 心理 学 研究 的 一 大 热点 ， 但 许多 密集 追踪 的 研究 分 析 仍 停留 在 较为 传统 的 方法 。 方 法 学 领域 已 涌现 出 
较 多 用 于 密集 追踪 数据 分 析 的 模型 方法 ， 较 为 主流 的 模型 包括 以 动态 结构 方程 模型 (Dynamic Structural 
Equation Model，DSEMD) 为 代表 的 自 上 而 下 的 建 模 方法 ， 以 及 以 组 和 迭代 多 模型 估计 (Group Iterative Multiple 
Model Estimation, GIMME) 为 代表 的 自 下 而 上 的 建 模 方法 。 二 者 均 可 以 方便 地 对 密集 追踪 数据 中 的 自 回 归 及 交 
又 滞后 效应 进行 建 模 。 

关键 词 ”密集 追踪 ， 时 间 序 列 ， 动 态 结 构 方 程 模型 ， 组 和 迭代 多 模型 估计 
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1 引言 常常 被 用 来 对 个 体 的 行为 以 及 状态 随时 间 变 化 的 
动态 过 程 进行 更 精细 的 测量 , 测量 内 容 包括 情绪 、 
治疗 或 自然 发 后 的 事件 中 个 体 的 相 
€ ss im : ie 人 格 、 认 知 和 行为 模式 等 方面 (Hofmans et al., 2019; 
" " i Ti ye "- is 的 p dk dE Lazarevic et al., 2020; Trull & Ebner-Priemer, 2014). 

* di KM tendi wi " 根据 Bringmann 等 人 (2017) 的 定义 , 密集 追踪 数据 
是 指 研 究 者 在 较 短 的 时 间 间 陋 内 ， 收 集 到 的 个 休 
FIT TA B6 — FE TE. D F $ Xt TAE zi St fz "BOX et as SEES 
行 描述 或 解释 ， 但 难以 揭示 短 时 间 内 个 体内 部 心 i oe RE 
理 状态 的 动态 变化 , 也 不 能 详细 解释 事件 对 个 体 m ges ar pent n 
KR de iri : ud Ted 息 ， 实 时 地 记录 被 试 随 时 间 变 化 的 动态 过 程 
d ie aus iid ue is s (Bolger & Laurenceau, 2013), 进一步 解释 个 体内 部 
2019)， 可 能 导致 开 究 者 对 干预 效果 及 其 潜在 机 制 心理 状态 和 行为 的 动态 变化 等 ( 张 银 普 等 ， 2016). 


产生 错误 解读 ( 张 银 普 S5, 2016; Setodji et al., 近年 来 ， 随 着 手机 等 移动 通讯 设备 的 发 展 和 
2019)。 国 此 ,为 了 更 深入 地 从 个 体内 的 角度 研究 。 使 用 ,研究 者 可 以 更 加 方便 、 频 繁 地 收集 有 关 被 


变量 间 的 关系 ,研究 者 往往 需要 采用 密集 扎 踪 方 。 试行 为 、 所 处 环境 、 心 理 和 生理 状态 等 方面 的 数 


法 , 更 密集 地 测量 个 体 状态 发 展 和 变化 的 细节 据 ; 数据 收集 过 程 的 便捷 化 以 及 收集 场景 的 多 样 
Eoo RAS E 化 等 ,也 都 有 利于 研究 者 探究 关于 不 同 场景 和 不 
密集 追踪 数据 有 测量 问 隔 得、 频次 高 的 特点 ， 同时 间 点 的 问题 。 如 今 密集 追踪 数据 已 经 被 广泛 

应 用 到 心理 学 的 多 个 领域 (Chen & Zhang, 2020; 

收 稿 日 期 ;2020.08.25 Windt e al., 2018), An, WRTA (Lanza et al., 
* 国家 自然 科学 基金 项 目 (31871128); 教育 部 人 文 社会 2014; Shiffman, 2009; Weinstock et al., 2017)、 压 力 


科学 研究 规划 基金 项 目 (18YJA190013)。 症状 (Chun, 2016; van de Schoot et al., 2015) 、 心 理 
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2014)、 饮 食 障碍 问题 (Munsch et al., 2009), 、 亲 密 
关系 (Howland & Rafaeli, 2010)、 人 格 障 碍 (Trull et 
al, 2008; Wright et al.，2015，2016) 和 语言 使 用 
(Lazarevic et al., 2020)55 , 

Trull 和 Ebner-Priemer (2013) 指 出 ， 密 集 追 踪 
方法 解决 了 很 多 回顾 性 调查 和 实验 室 方法 的 局 限 ， 
如 潜在 的 记忆 偏差 、 对 个 体 日 常生 活 场 景 的 还 原 
度 不 足 等 (Bolger & Laurenceau, 2013; Sened et al., 
2018)。 考 虑 到 密集 追踪 方法 存在 的 诸多 优势 ， 国 
内 一 些 研究 者 开始 使 用 密集 追踪 数据 来 验证 实验 
假设 ， 如 探究 干预 方法 对 个 体 心 理 和 行为 变化 的 
影响 ( 陈 明 瑞 ， 周 萍 ，2017; 郑 文 倩 等 , 2019)、 环 
境 因 素 对 个 体 身 心 状态 的 影响 ( 安 媛 媛 ， 徐 奈 ， 
2017)、 个 体 在 日 常生 活 中 的 压力 应 对 ( 邵 华 等 ， 
2019) 以 及 个 体 在 休闲 活动 中 主观 幸福 感 的 变化 
等 ( 任 杰 等 , 2010; 张 峰 ,， 耿 晓 伟 ,，2015)。 这 些 研 
究 多 采用 较为 传统 的 多 层 线性 模型 (Multilevel 
Model, MLM) 对 密集 追踪 数据 进行 分 析 ( 安 姐 姐 ， 
IREL, 2017; 邵 华 等 ,2019; 张 银 普 等 , 2016)。 然 
而 ,密集 追踪 方法 的 优势 更 多 地 在 于 探究 变量 在 
变化 过 程 中 的 动态 效应 ( 即 变量 随时 间 的 变化 可 
能 受到 变量 本 身 或 其 他 外 部 因素 的 影响 ; Hamaker 
& Wichers, 2017; Molenaar, 2004; Molenaar & 
Campbell，2009)， 因 此 对 密集 追踪 数据 的 分 析 方 
法 提出 了 更 高 的 要 求 。 

为 了 帮助 研究 者 们 增进 对 于 密集 追踪 数据 方 
法 的 了 解 ， 下 文 将 首先 对 常见 的 密集 追踪 研究 方 
法 及 其 优势 进行 整体 梳理 ， 然 后 介绍 密集 追踪 数 
据 中 较为 常用 的 分 析 方 法 (MLM)， 以 及 目前 较 新 
的 分 析 方 法 ,主要 包括 动态 结构 方程 模型 
(Dynamic Structural Equation Model, DSEM) 和 组 
迭代 多 模型 估计 (Group Iterative Multiple Model 
Estimation, GIMME)。 为 了 更 直观 地 比较 MLM, 
DSEM 和 GIMME 在 实际 建 模 过 程 中 的 特点 ， 本 
文 分 别 利用 这 三 种 模型 方法 进行 实证 的 演示 和 比 
较 。 最 后 , 为 了 帮助 研究 者 针对 自身 研究 的 实际 
情况 灵活 选择 适合 的 模型 方法 ， 本 文 还 将 对 各 种 
主流 模型 方法 的 优 缺 点 和 适用 范围 进行 总 结 ， 以 
供 参 考 。 


2 密集 追踪 研究 方法 及 优势 


2.1 密集 追踪 研究 方法 
密集 追踪 方法 通常 又 称 动态 评估 法 


uu 


(Ambulatory Assessment, AA; Bolger et al., 2003; 
Sened et al., 2018)， 其 另外 几 种 常见 的 术语 有 : (1) 
经 验 取样 法 (Experience Sampling Method, ESM; 
Csikszentmihalyi & Larson, 2014), 广泛 应 用 于 人 
格 和 社会 心理 学 领域 ， 主 要 使 用 主观 问卷 进行 测 
量 ; (2) 日 记 法 (Daily Diary; Bolger & Laurenceau, 
2013) 和 日 常 过 程 研究 法 (Daily Process Approach; 
Tennen et al.,， 2000), 通过 被 试 的 每 日 记录 ,追踪 
有 关 变 量 随时 间 变 化 的 过 程 ; (3) 生 态 瞬 时 评估 法 
(Ecological Momentary Assessment, EMA; Stone & 
Shiffman, 1994), 常用 于 临床 和 药物 研究 。 典 型 的 
EMA 研究 收集 被 试 在 日 常生 活 中 每 日 特定 时 刻 
的 状态 和 行为 数据 ， 在 这 些 特定 的 应 答 时 刻 中 ， 
被 试 或 是 被 动 接受 来 自 外 界 的 提示 , 或 是 自主 发 
起 某 些 行为 (如 ， 被 试 需要 在 感受 到 对 成 瘾 药物 的 
强烈 欲望 后 进行 作答 ; Lazarevic et al., 2020); (4) 
动态 监测 法 (Ambulatory Monitoring; e.g., O'Brien 
et al.，2000), 被 试 需要 佩戴 特定 的 便携 式 仪器 ， 
对 特定 的 检测 指标 (如 血压 ) 进 行 密集 测量 。 

一 般 来 说 ， 密 集 追 踪 研 究 方法 可 以 被 广义 地 
分 为 三 类 (Fisher & March, 2012; Trull & Ebner- 
Priemer，2013), 包括 自我 报告 式 的 自我 评估 ; 观 
察 式 的 自我 评估 ; 生理 、 心 理 或 行为 的 自我 评估 。 
在 自我 报告 式 的 评估 中 ,根据 不 同 的 实验 目的 ， 
被 试 可 以 在 固定 /随机 间隔 的 时 间 点 进行 自我 评 
fh, 或 者 在 特定 事件 发 生 后 进行 作答 (Reis et al., 
2014)。 不 同 的 自我 报告 方式 各 有 优势 和 不 足 : (1) 
国定 时 间 间 隔 的 自我 报告 方式 要 求 被 试 在 约定 好 
的 否 干 时 间 点 回答 相应 的 问题 (Fisher & March, 
2012)， 其 数据 收集 的 过 程 不 会 突然 打扰 到 被 试 ， 
数据 收集 方式 本 身 简单 易 操作 。 然 而 ,固定 时 间 
间隔 的 自我 报告 也 可 能 因为 只 考虑 固定 时 间 点 发 
生 的 事件 ， 而 忽略 了 其 他 时 间 点 的 事件 可 能 产生 
的 影响 ; 同时 ， 当 测量 时 间 间 隔 过 大 ,也 可 能 使 
得 被 试 出 现 记忆 衰退 的 现象 (Fisher & March, 
2012; Hofmans et al., 2019). (2) 在 随机 时 间 间 隔 的 
自我 报告 中 ， 被 试 需要 不 定期 地 配合 研究 者 完成 
数据 的 收集 (Fisher & March，2012)， 因 此 可 以 随 
机 抽 到 一 些 具 有 代表 性 的 经 历 ; 但 是 对 于 被 试 而 
言 ， 每 次 数据 收集 的 时 间 是 未 知 的 ， 这 种 不 确定 
性 可 能 对 被 试 的 日 常生 活 带 来 较 大 的 影响 和 负担 
(Fisher & March, 2012; Hofmans et al., 2019). (3) 
关注 特定 事件 的 数据 收集 方式 ， 是 指 被 试 在 发 生 
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特定 事件 后 主动 配合 研究 者 完成 数据 的 收集 
(Fisher & March，2012)， 其 优势 在 于 可 以 收集 到 
与 研究 者 感 兴 趣 的 事件 相关 的 数据 。 然 而 ， 这 种 
方法 收集 到 的 数据 取决 于 被 试 对 事件 的 反应 ， 因 
此 研究 者 需要 事先 对 被 试 进行 训练 , 使 其 能 意识 
到 特定 事件 的 发 生 并 主动 及 时 向 研究 者 报告 
(Fisher & March, 2012; Hofmans et al., 2019)。 

另外 ,与 自我 报告 的 评估 不 同 ， 观 察 性 自我 
评估 和 生物 、 心 理 或 行为 的 自我 评估 主要 收集 客 
观 数据 和 指标 (如 ,被 试 所 处 环境 中 的 声音 、 被 试 
运动 的 速度 ， 以 及 心率 等 生理 方面 的 指标 数据 )， 
而 不 是 主观 报告 的 数据 (Trull & Ebner-Priemer, 
2013)， 因 此 很 少 使 被 试 产生 心理 负担 ， 同 时 其 数 
据 收 集 也 不 易 受 到 被 试 的 主观 影响 。 具 体 关 于 不 
同 密集 追踪 方法 的 使 用 可 以 参考 Fisher 和 March 
(2012)、van Roekel 等 人 (2019) 的 文章 ， 以 及 本 文 
关于 密集 追踪 研究 中 的 实际 研究 示例 。 
2.2 ”密集 追踪 研究 方法 的 优势 

相 比 于 传统 横断 研究 的 方法 ， 密 集 追 踪 研 究 
方法 的 优势 主要 体现 在 : (]) 与 横断 研究 中 回顾 式 
自我 报告 的 方式 相 比 ， 密 集 追 踪 方 法 对 变量 的 测 
量 更 及 时 、 更 接近 被 试 的 真实 经 历 ， 可 以 更 好 地 
减少 测量 过 程 中 出 现 的 记忆 偏差 (Stone et al, 
1999)。 因 此 ， 所 产生 的 密集 追踪 数据 对 于 事件 的 
描述 往往 更 加 准确 ,更 能 代表 被 试 的 真实 反应 
(Sened et al., 2018)。(2) 与 横断 研究 中 的 实验 室 方 
法 相 比 ， 密 集 追 踪 方 法 提供 了 多 个 时 间 点 和 场景 
的 测量 , 使 得 研究 者 可 以 探究 个 体内 部 心理 构 念 
随时 间 的 动态 发 展 及 相互 作用 ( 张 银 普 等 ，2016; 
Trull & Ebner-Priemer, 2014)。(3) 密 集 追 踪 方 法 的 
测量 随 着 被 试 的 日 常生 活 情境 进行 ,不仅 局 限于 
实验 室 环境 , 极 大 降低 了 社会 期 望 以 及 被 试 自 我 
暴露 意愿 对 实验 数据 的 影响 ,因此 有 更 好 的 生态 
效 度 和 外 部 效 度 (Sened et al., 2018). (4) 5 4218 E: 
方法 对 多 变量 重复 测量 ， 人 允许 研究 者 分 别 探究 个 
体 间 和 个 体内 的 变异 (Curran & Bauer, 2011; 
Sened et al., 2018)。 许 多 心理 学 研究 的 理论 和 假设 
主要 针对 个 体内 部 的 过 程 (Intraindividual 
Processes), 但 其 数据 分 析 的 本 质 往往 是 对 个 体 间 
关系 (Interindividual Processes) 的 探讨 (如 ，Lau 等 
人 (2016) 针 对 个 体 间 自 变量 水 平 变 化 对 因 变 量 的 
影响 ， 得 出 “孕妇 对 压力 的 积极 感知 会 促进 其 积 
极 应 对 的 策略 ， 进 而 提高 身心 健康 ”这 一 个 体内 


的 结论 )， 因 此 这 样 得 到 的 关于 个 体内 关系 的 结论 
可 能 是 存在 问题 的 (Curran & Bauer, 2011). 

相 比 于 传统 的 追踪 方法 ， 密 集 追 踪 也 具有 其 
独特 的 优势 。 传 统 追 踪 方 法 收集 的 数据 重复 观察 
的 数量 较 少 (通常 少 于 10 次 ; Jahng, 2008)， 并 且 测 
量 时 间 间 隔 较 大 ， 因 此 主 要 用 来 探究 在 相对 宏观 
的 时 间 尺 度 上 ， 因 变量 总 体 均值 随时 间 的 变化 
(McNeish & Hamaker, 2019; Nesselroade, 1991; 
Setodji et al., 2019), 但 不 适合 用 于 研究 个 体内 变 
量 的 动态 变化 过 程 (Jahng, 2008). 而 密集 追踪 方法 
对 变量 的 测量 间隔 和 持续 时 间 更 短 、 频 次 更 多 
(Hamaker & Wichers, 2017; Walls & Schafer, 2006), 
可 以 更 深入 地 反映 变量 随时 间 变 化 的 复杂 模式 
(Bolger & Laurenceau, 2013)， 人 允许 研究 者 对 被 试 
的 行为 或 心理 社会 过 程 随时 间 产 生 的 细微 变化 ， 
以 及 相关 的 共 变 机 制 等 进行 研究 (Collins, 2006)。 
相 比 于 传统 的 追踪 方法 ， 密 集 追 踪 方 法 在 被 试 内 
水 平 的 数据 点 较 多 ， 有 足够 的 检验 力 检验 被 试 内 
水 平 的 研究 假设 , 改善 了 传统 追踪 方法 中 被 试 内 
水 平 的 数据 点 较 少 而 统计 检验 力 不 足 的 问题 
(Fraley & Hudson, 2014; Hofmans et al., 2019)。 

总 的 来 说 ， 密 集 追 踪 研 究 方法 能 更 好 地 保证 
数据 的 有 效 性 和 真实 性 (Sened et ab, 2018), 允许 
研究 者 从 个 体 的 角度 探究 变量 的 变化 关系 , 深入 
分 析 个 体内 部 的 心理 变化 过 程 (Baltes & 
Nesselroade, 1979; Molenaar, 2004)。 利 用 密集 追踪 
方法 , 研究 者 可 以 建立 更 有 针对 性 的 、 针 对 于 个 
体 的 模型 ， 从 而 更 加 精准 地 识别 个 体 特征 并 对 其 
进行 更 有 效 的 治疗 (Wright & Woods, 2020); 并 通 
过 探究 变量 在 被 试 间 水 平 的 差异 , 检验 单个 个 体 
的 效应 是 否 能 够 推广 或 者 泛 化 到 更 多 的 样本 中 
(Hofmans et al., 2019)。 


3 密集 追踪 数据 的 传统 分 析 方法 


密集 追踪 数据 具有 上 典型 的 分 层 结构 ， 即 水 平 
1( 时 间 点 水 平 ) 的 基本 单元 嵌 套 于 水 平 2( 个 体 水 平 ) 
的 结构 ,水 平 2 的 结构 可 以 继续 般 套 于 更 高 的 水 
平 (如 ,组 别 水 平 )， 因此 传统 上 利用 MLM 进行 分 
析 ( 张 银 普 45, 2016; Schwartz & Stone, 2007; 
Walls & Schafer, 2006), 将 变量 的 变异 来 源 分 为 个 体 
内 和 个 体 间 等 多 个 水 平 (Walls & Schafer, 2006)。 

3.1 方法 原理 
以 两 层 的 MLM 为 例 , 一 般 来 说 ,密集 追踪 数 
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据 中 的 水 平 1 关注 个 体内 变量 在 不 同时 间 点 之 间 
的 关系 ， 而 水 平 2 关注 个 体 水 平 上 的 变量 间 的 关 
系 。 对 于 水 平 1( 时 间 点 水 平 ), 以 存在 一 个 水 平 1 
预测 变量 的 情况 为 例 , 个 体 i 在 测量 时 间 点 t 的 
模型 公式 如 下 (Heck & Thomas, 2015): 

Y, = Poi + BiXi + En (1) 
HP, X, 表示 个 体 i 在 第 t 个 测量 时 间 点 的 测量 ， 


(Hamaker & Grasman, 2015): 
®© = © +u; (5) 
其 中 ，@ 表示 所 有 个 体 平均 的 自 回归 效应 ，w; 表 
示 个 体 i 的 自 回归 效应 与 平均 自 回 归 效 应 之 间 的 
m, RAES, u; ~N (0, 02, )。 
3.2 ”传统 多 层 线 性 模型 的 不 足 
传统 的 MLM 可 以 解释 不 同 个 体 间 变量 随时 


6 和 分 别 表示 截 距 和 和 斜率 ，si 表示 误差 项 
(er~ N(0,07)). 另外 ,在 密集 追踪 数据 的 模型 中 ， 
时 间 太 ,可 以 作为 预测 变量 加 入 到 水 平 1 的 模型 中 。 
对 于 水 平 2( 个 体 水 平 ), 同样 以 存在 一 个 水 平 

2 预测 变量 的 情况 为 例 ， 其 模型 表达 式 如 下 : 
Boi = Yoo + Noi + to; (2) 
Bu = Yor + Wi + (3) 
Rh, fy; 和 Bi; 分 别 为 水 平 1 中 的 随机 截 距 和 随 
机 和 斜率， 到 为 水 平 2 的 预测 变量 ，yi。、i 分 别 
为 水 平 2 的 预测 变量 万 对 Bo, 、 忆 的 斜率 ，7o0 、 
Yo 表示 对 应 的 截 距 。 各 变量 在 水 平 1 的 斜率 和 截 
距 可 能 在 不 同 个 体 间 存在 变异 (随机 效应 )， 且 这 
一 变异 可 能 可 以 被 水 平 2 的 预测 变量 解释 。 另 外 ， 
水 平 1 和 水 平 2 的 模型 均 可 以 处 理 含 有 多 个 预测 


间 变 化 的 关系 曲线 存在 的 差异 (Walls & Schafer, 
2006), 但 是 这 种 模型 不 够 灵活 。 在 分 析 动 态 变 化 
效应 、 探 究 变 量 间 关 系 的 强度 和 方向 随时 间 变 化 
的 趋势 时 ,传统 的 MLM 仍 存在 一 定 的 局 限 (Li 
et al., 2006)。 首先 ， 当 研究 问题 关注 的 是 个 体内 心 
理 过 程 的 动态 变化 ， 如 变异 性 和 不 稳定 性 时 ， 传 
统 的 MLM 方法 在 分 析 个 体内 发 展 模式 中 的 个 体 
间 差 异 时 可 能 会 产生 不 正确 的 参数 估计 (Jahng， 
2008)， 且 不 能 解决 变量 随时 间 变 化 的 轨迹 在 个 体 
间 存 在 较 大 差异 的 问题 (Shiyko et al., 2012)。 其 次 ， 
MLM 通常 假设 随机 效应 服从 一 定 的 分 布 类 型 ， 
如 正 态 分 布 ; 但 通常 情况 下 ,这 种 假设 很 难得 到 
满足 ， 因 此 基于 这 些 假设 得 到 的 模型 参数 也 可 能 
是 不 准确 的 (Piccirillo & Rodebaugh, 2019)。 针 对 


变量 的 情况 。 根 据 实际 研究 中 密集 追踪 数据 的 具 
体 情况 ,可 以 对 其 进行 更 多 层 的 建 模 ， 如 测量 时 
kB TMA, 个 体 进一步 舰 套 于 所 在 的 群 
体 等 。 

在 使 用 MLM 分 析 传 统 的 追踪 数据 时 , 通常 
假设 不 同 个 体 在 水 平 1 的 残 差 方差 o 均 相 等 , 但 
在 密集 追踪 数据 中 ， 这 种 假设 很 可 能 会 带 来 一 些 
问题 (Jahng & Wood, 2017)， 如 参数 标准 误 的 有 偏 
估计 等 (Hamaker & Grasman，2015)。 密 集 追 踪 数 
据 中 多 次 重复 测量 的 时 间 间 隔 较 短 , 个 体 在 当前 
时 间 点 的 测量 很 可 能 会 受到 先前 测量 的 影响 ， 因 
此 通常 需要 假设 水 平 1 的 残 差 带 有 自 回归 的 效应 
(Hamaker & Grasman, 2015)。 以 考虑 一 阶 自 回归 
效应 (Lag-1 Autoregressive Effects) 为 例 ， 此 时 式 
() 中 的 残 差 e, 可 以 表示 为 (Hamaker & Grasman, 
2015): 


I 


En = Q8; + en (4) 
其 中 ， 残 差 项 e, 可 以 视 为 模型 中 没有 考虑 到 的 因 
素 的 影响 ， 服 从 正 态 分 布 , 即 e -N (0,02). à 


这 一 问题 ,Ke 等 人 (2019) 指 出 , 将 贝 叶 斯 多 层 结 
构 方 程 模型 (Bayesian Multilevel Structural Equation 
Modeling) 与 更 稳健 的 、 使 用 非 正 态 分 布 ( 如 
Student's t 分 布 , Student's t Distribution; 混合 正 态 
分 布 ，Mixture of Normal Distributions; Tong & 
Zhang, 2012; Zhang, 2013; Zhang et al., 2013) 的 贝 
叶 斯 估计 方法 相 结合 ,可 能 可 以 解决 随机 效应 不 
符合 正 态 分 布 时 参数 估计 不 准确 的 问题 。 

另外 , 在 MLM 中 , 研究 者 往往 需要 通过 中 心 
化 来 分 解 组 间 水 平和 组 内 水 平 的 效应 (Curran & 
Bauer, 2011)。 传 统 上 对 于 预测 变量 的 处 理 方法 包 
括 个 体 均值 中 心 化 、 总 均值 中 心 化 等 , 但 这 些 方 
法 的 适用 范围 较为 有 限 ， 在 一 些 情况 下 容易 产生 
较 大 的 偏差 (Asparouhov & Muthén, 2019a; Curran 
& Bauer, 2011)。 如 , 在 存在 情境 效应 ! (Contextual 
Effects; Liidtke et al., 2008) 或 自 回归 效应 (Nickell, 
1981)， 以 及 预测 变量 为 (两 个 及 以 上 类 别 的 ) 分 类 
变量 且 存 在 非 线性 连接 函数 (Link Function) 的 情 


示 个 体 i 的 自 回归 效应 ， 即 该 个 体 在 当前 测量 时 
间 点 的 测量 偏差 ; 在 多 大 程度 上 受到 上 一 测量 时 
间 点 的 影响 。 自 回归 系数 @ 在 水 平 2 的 表达 式 为 


F 


Bl, 在 控制 了 水 平 1 的 自 变量 x, 的 情况 下 , 水 平 2 的 平 
均 水 平 X, 仍 对 因 变 量 存在 额外 的 效应 。 
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况 下 ,传统 方法 均 会 产生 较 大 的 参数 估计 偏差 
(Asparouhov & Muthén, 2019a)。 当 缺失 数据 是 随 
机 缺失 (Missing at Random，MAR) 而 不 是 完全 随 
机 缺失 (Missing Completely at Random, MCAR) 时 ， 
MLM 中 的 样本 均值 无 法 体现 预测 变量 真实 的 均 
值 ， 传 统 的 中 心 化 处 理 同样 可 能 带 来 参数 估计 上 
的 偏差 (Asparouhov & Muthén, 2019a)。 


4 密集 追踪 数据 分 析 方 法 的 新 进展 


为 了 克服 传统 方法 在 密集 追踪 数据 中 的 局 限 
性 ,帮助 研究 者 更 好 地 挖掘 个 体内 部 动态 变化 的 
过 程 和 作用 机 制 ， 下 面 将 主要 介绍 两 种 较 新 的 密 
集 追 踪 数 据 的 分 析 方 法 ,包括 DSEM (Asparouhov 
et al.，2018) 和 结合 了 个 体 和 群体 信息 的 GIMME 
(Gates & Molenaar，2012)。 根 据 Hamaker 等 人 
(2018) 的 界定 ， 多 个 体 的 密集 追踪 数据 分 析 方法 
可 以 分 为 自 上 而 下 的 方法 (Top-Down Approach) 
和 自 下 而 上 的 方法 (Bottom-Up Approach) WK 


4.1.1 方法 原理 

DSEM 可 以 认为 是 动态 因子 模型 (Dynamic 
Factor Model, DFM; Ferrer & Nesselroade, 2003; 
Molenaar, 1985; Molenaar et al., 1999; Zhang & 
Nesselroade，2007) 在 多 个 个 体 情 况 下 的 扩展 
(Asparouhov et al., 2018)。 它 结合 多 层 模 型 的 思路 ， 
将 第 i 个 个 体 在 第 t 个 时 间 节 点 的 观测 变量 多 ， 
分 解 为 个 体 i 的 特异 贡献 二;,、 时 间 t 的 特异 贡献 
Y,,, 以 及 除去 个 体 和 时 间 的 特异 贡献 后 , 个 体 i 
在 时 间 节 点 f 的 波动 变化 7; 这 三 个 成 分 ， 即 
Y, =Y , + Y; ; +Y; , (Asparouhov et al., 2018)。 例 如 ， 
对 被 试 群体 在 工作 日 中 感受 到 的 工作 压力 和 积极 
情绪 水 平 进行 测量 ， 其 中 个 体 ?z 在 时 间 节 点 上 的 工 
TER FACE Y, , 将 同时 取决 于 个 体 工 在 工作 压力 
上 的 一 般 基线 水 平 交 .时 间 节 点 ! 对 个 体 工作 压 
力 的 一 般 效 应 Y, (如 ,一 般 来 说 ,每 到 周一 时 个 
体 的 工作 压力 都 会 增高 )， 以 及 除去 个 体 和 时 间 的 
一 般 性 效应 后 , 个 体 i 在 时 间 节 点 1 的 工作 压 


类 。 自 上 而 下 的 方法 通过 较 高 水 平 的 模型 向 下 对 
较 低 水 平 的 动态 效应 进行 解释 (如 MLM，DSEM 
等 )， 这 里 的 高 低 次 序 均 为 相对 的 而 不 是 绝对 的 
高 低 ( 如 在 三 水 平 的 MLM 中 ,水平 2 相对 于 水 平 
1 是 较 高 水 平 , 但 相对 水 平 3 是 较 低 的 水 平 )。 而 


AE 4E Y, a (如 公式 (6)~(7) 所 示 ， 这 些 变化 的 来 源 可 
能 包括 : 该 个 体 在 该 时 间 节 点 的 积极 情绪 变化 ， 
先前 时 间 点 的 积极 情绪 和 工作 压力 水 平 对 当前 时 
间 点 的 影响 等 )。 

最 完整 的 DSEM i HEAR LTR Yo a Y, Wh 


自 下 而 上 的 方法 首先 对 每 个 个 体 的 数据 进行 分 
析 ， 再 进一步 向 上 搜寻 不 同 个 体 间 相 似 的 模式 
(如 GIMME), 
41 自 上 而 下 的 方法 : 动态 结构 方程 模型 

DSEM 承接 了 传统 MLM 自 上 而 下 的 思想 ， 
同时 ， 由 于 它 采 用 了 贝 叶 斯 的 框架 ， 因 此 相 比 传 
Zi; MLM 更 为 灵活 (Asparouhov et al., 2018).DSEM 
结合 了 多 种 模型 方法 的 优势 , 包括 MLM, 、 时 间 序 
列 模型 (Time-Series Modeling), 、 结 构 方 程 模型 
(Structural Equation Modeling，SEM) 以 及 时 变 效 
应 模型 (Time-Varying Effects Modeling, TVEM), 
被 认为 是 用 于 分 析 密 集 追 踪 数 据 最 有 前 景 的 方法 


BOY, 这 三 个 成 分 ， 但 不 是 所 有 的 密集 追踪 研究 
均 存 在 时 间 特 异 的 效应 ， 因 此 不 是 所 有 的 密集 追 
踪 人 研究 均 适合 采用 最 完整 的 DSEM。 如 , 在 密集 
追踪 的 研究 中 ， 如 果 1=1 MRA PAW BPM 
测 时 间 点 的 测量 ,而 不 是 某 个 特定 的 时 间 节 点 ( 比 
如 t=1 都 表示 每 周 的 第 一 天 ) 时 ，t=1 对 于 每 个 个 
体 来 说 没有 一 致 的 含义 ， 因 此 数据 本 身 很 可 能 不 
存在 时 间 特 异 的 效应 (Asparouhov et al., 2018). 
Asparouhov 等 人 (2018) 指 出 ， 研 究 者 可 以 根据 实 
际 研究 的 需要 而 舍弃 马 , 这 一 成 分 只 考虑 个 体 
水 平 的 随机 效应 ， 即 采用 两 层 的 DSEM (Two- 
Level DSEM); 甚至 进一步 舍弃 及,， 只 考虑 单 


之 一 (Zhou et al., 2020)。 它 强调 变量 在 时 间 尺 度 上 
的 动态 变化 过 程 ， 这 里 的 “动态 "可 以 理解 为 变量 
在 当前 时 刻 的 测量 会 受到 变量 自身 在 先前 时 间 点 
的 影响 ( 自 回归 效应 )， 以 及 受到 其 他 变量 在 先前 
(交叉 滞后 效应 ，Cross-lagged Effects) 以 及 当前 时 
间 点 (同时 效应 ，Contemporaneous Effects) 的 测量 


的 影响 。 另 外 , 在 DSEM 中 , 各 变量 的 变异 均 可 
通过 个 体 特异 和 时 间 特 异 的 变量 进行 解释 。 


层 的 DSEM 模型 (Single-Level DSEM)。 两 层 和 
单 层 的 DSEM 均 为 总 框架 下 的 特殊 的 情况 ， 
此 ,本文 将 对 总 的 DSEM 框架 下 的 方法 原理 进行 


介绍 


首先 , 与 DFM XW, DSEM 在 水 平 1 考虑 了 
变量 在 时 间 节 点 t 的 测量 可 能 受到 当前 及 先前 工 
次 测量 的 影响 , 包括 测量 模型 和 结构 方程 模型 
(Asparouhov et al., 2018): 
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L L 
F, =y + > Ay th it-1 + > Rolin t 
1-0 1-0 


L 


> Ky Xiii + Ei (6) 
1-0 
L L 
Mh, it =%, in + > By uh it-1 + ` QuY, i, -1 + 
1-0 1-0 
L 
> Vy it X 1, 41-1 + Gi (7) 
1-0 


RB, Ma RIR Y a SOTA, A 为 因子 负 
BEM, R, K., B. 、O 、 工 表示 回归 系数 和 矩阵 。 
Eins ii 表示 残 差 ,均值 均 为 0。 另外 , DSEM HE 
架 在 水 平 2 BACK 3 的 模型 公式 如 下 (Asparouhov 
et al., 2018): 


Y, i =v +A, i +K Xa ite; (8) 
7,;705*Byp i TX i téa i (9) 
Y, =v; t Aq + KyXS +E, (10) 
Th, = 03 By, *V4X4, 55, (11) 


其 中 , 及; 为 一 组 个 体 特异 (在 不 同 个 体 间 存 
在 差异 , 但 不 随时 间 变 化 ) 的 变量 , X, 为 一 组 时 
间 特 异 (在 不 同时 间 节 点 间 存 在 差异 , 但 不 随 个 
体 变化 ) 的 变量 。72, CFI ms LAT BIO Y, UY, E 
的 潜 变 量 。 总 的 来 说 , DSEM 在 水 平 1 考虑 了 时 间 
节点 上 的 测量 与 当前 及 先前 工 次 测量 的 回归 关系 ， 
EKE 2 考虑 个 体 特异 的 效应 ; 而 且 ， 与 传统 多 
层 模型 不 同 的 是 , DSEM 还 在 水 平 3 考 虑 了 时 间 特 
异 的 效应 。 

具体 来 说 , DSEM 在 水 平 1 中 的 回归 系数 、 因 
子 负载 和 截 距 项 等 参数 的 随机 效应 可 以 进一步 在 
水 平 2 和 水 平 3 进行 解释 ,同时 考虑 不 同 个 体 和 
时 间 节 点 的 特异 贡献 (Asparouhov et al., 2018)。 对 
FKE 1 的 系数 s (如 ，s 可 以 指 代 回 归 系 数 和 矩阵 
Ry, 中 的 元 素 、 因 子 负载 矩阵 Au 六 中 的 元 素 、 截 
距 项 a, 等 ) 其 随机 效应 可 以 在 水 平 2 和 水 平 3 
被 分 解 为 s,; 和 95， 两 个 成 分 ， 并 分 别 被 个 体 特 
异 和 时 间 特 异 的 变量 解释 (Asparouhov_ et al., 
2018)。 如 , 在 最 完整 的 DSEM 框架 下 , 个 体积 极 
情绪 的 自 回归 系 数 的 随机 效应 将 被 分 解 为 个 体 和 
时 间 两 个 方面 的 来 源 。 这 种 随机 效应 在 水 平 2 可 
能 受到 人 格 特质 (个 体 特异 的 变量 ) 的 影响 ， 同 时 
在 水 平 3 受到 特定 工作 日 (时 间 特 异 的 变量 ) 的 影响 。 
类 似 地 ， 对 于 模型 中 的 残 差 ， 其 方差 的 随机 效应 也 可 
以 抽象 地 表示 为 指数 函数 的 形式 VV = Exp(s ,+s ,). 


从 而 可 以 分 别 利 用 个 体 和 时 间 特 异 的 变量 对 其 随 
机 效应 进行 解释 (Asparouhov et al., 2018)。 

McNeish 和 Hamaker (2019) 指 出 ， 相 比 于 
MLM, DSEM 的 亮点 在 于 它 可 以 很 方便 地 考虑 变 
量 的 自 回归 效应 和 交叉 滞后 的 效应 。 标 准 的 MLM 
软件 也 可 以 创造 滞后 时 间 点 的 新 变量 加 入 模型 ， 
从 而 进行 和 DSEM 类 似 的 分 析 ; 但 当时 间 间 隔 不 
等 或 存在 缺失 数据 时 , MLM 中 的 缺失 数据 点 也 随 
之 增多 ， 从 而 对 参数 估计 带 来 更 严重 的 问题 
(McNeish & Hamaker, 2019)。 而 DSEM 是 利用 卡 
尔 曼 滤波 (Kalman Filter; Harvey, 1989), 基于 上 一 
时 间 点 的 观测 对 当前 缺失 时 间 点 的 数据 迭代 地 进 
行 预测 和 更 新 ,不 需要 依赖 按 列 删除 (Listwise 
Deletion)， 因 此 在 处 理 缺失 数据 方面 有 较 大 的 优 
HA, Ah, WN 3.2 节 所 述 , 在 MLM F, 人 研究 者 往 
往 需要 通过 中 心 化 来 分 解 组 间 水 平和 组 内 水 平 的 
效应 (Curran & Bauer，2011), 但 传统 的 显 变 量 中 
心 化 在 密集 追踪 数据 中 表现 不 佳 (Asparouhov & 
Muthén, 20193), DSEM 允许 对 潜 变 量 进行 建 模 ， 
采用 潜 变 量 中 心 化 的 方法 ( 即 在 结构 方程 模型 前 
at, HEE X BAA Xia Xar X4 T A AY 
别 进行 建 模 ) 来 排除 来 源 于 个 体 和 时 间 的 一 般 性 
影响 (Asparouhov et al., 2018)。 在 大 多 数 情 况 下 ， 
这 种 潜 变 量 中 心 化 的 方法 在 参数 估计 方面 的 表现 
均 比 传统 显 变 量 中 心 化 的 方法 更 优 ， 能 很 好 地 解 
决 传统 的 显 变 量 中 心 化 带 来 的 问题 (Asparouhov 
& Muthén, 2019a)。 

需要 注意 的 是 , DSEM 要 求 时 间 节 点 上 为 整数 ， 
但 在 实际 数据 中 这 一 条 件 很 可 能 无 法 满足 
(Asparouhov et al., 2018)。 例 如 , 各 测量 节点 之 间 
的 时 间 间 隔 可 能 参差 不 齐 (West，2018), 使 得 时 间 
节点 1 不 为 整数 ,为 了 解决 这 一 问题 ,DSEM 需要 
建立 另 一 个 度量 网 络 ， 对 时 间 节 点 进行 转换 和 处 
理 ， 从 而 得 到 一 系列 近似 的 整数 值 (Asparouhov et 
al.，2018)。 例 如 ， 当 四 个 测量 节点 的 间隔 分 别 为 
0.5 天 、0.9 天 和 2.1 天 ,可 以 选择 一 个 较 小 的 时 间 
间隔 6 = 0.5 作为 网 格 的 最 小 单位 ， 将 四 个 测量 时 
间 节 点 近似 转化 为 0= 0, 1= 1,4=2,t4=4。 由 于 
不 存在 观测 值 , DSEM 将 该 节点 视 为 缺失 值 处 理 
(Asparouhov et al., 2018)。Asparouhov 等 人 (2018) 
指出 ， 过 小 的 时 间 间 隔 5 会 导致 较 多 的 缺失 值 ， 
而 间隔 偏 大 则 可 能 会 导致 多 个 时 间 节 点 转化 为 同 
一 个 整数 ， 因 此 需要 慎重 考虑 6 的 设置 。 对 于 如 
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何 选 择 较为 合适 的 时 间 间 隔 6 值 ，Asparouhov 等 
人 (2018) 提 供 了 较为 详细 的 建议 。 

一 般 来 说 ， 贝 叶 斯 估计 不 需要 依赖 渐进 理论 ， 
在 小 样本 的 情况 下 相 比 频率 学 派 的 方法 具有 潜在 
的 优势 (van de Schoot et al., 2014), 但 这 种 优势 不 
是 仅 通过 单纯 改变 估计 方法 便 能 得 以 实现 的 
(Depaoli & Clifton, 2015; van de Schoot et al., 
2015), 往往 还 需要 研究 者 结合 实际 情况 对 参数 设 
置 合适 的 先 验 分 布 。 与 其 他 许多 贝 叶 斯 框架 下 的 
模型 类 似 (Depaoli & van de Schoot, 2017; McNeish, 


化 , 还 可 以 考虑 变量 间 关 系 随时 间 的 动态 变化 
(Asparouhov et al., 2018; Zhou et al., 2020)。 另外 ， 
DSEM 对 模型 参数 的 估计 采用 贝 叶 斯 模型 的 框架 ， 
相 比 于 传统 的 频率 学 派 方法 ， 贝 叶 斯 方法 对 于 复 
杂 模 型 的 处 理 更 有 优势 ( 张 沥 今 等 ,2019; Zhou 
et al., 2020)。Song 和 Zhang (2014) 曾 将 多 层 结 构 
方程 模型 与 DFM 相 结合 ， 提 出 多 层 的 动态 因子 模 
型 (Multilevel Dynamic Factor Model, MDFM), 1H 
当 MDFM 中 存在 随机 效应 时 ,往往 会 导致 模型 无 
法 估计 的 问题 。 而 DSEM 对 于 测量 节点 和 随机 效 


2017; Zondervan-Zwijnenburg et al., 2017), 在 样 


应 的 设置 没有 数量 的 限制 ， 即 使 存在 较 多 的 测量 


本 量 较 小 时 使 用 Mplus 默认 的 无 信息 先 验 会 使 得 
DSEM 更 倾向 于 产生 有 偏 的 估计 结果 (McNeish, 
2019)。 因 此 , 在 样本 量 较 小 时 ,建议 研究 者 尽 可 
能 结合 前 人 研究 和 其 他 先 验 知识 对 参数 设置 有 信 
息 的 先 验 分 布 ; 在 缺乏 先 验 信 息 的 情 识 下 ， 即 使 
是 根据 参数 取 值 范围 的 直观 判断 来 对 参数 设置 容 
许 范 围 限 制 的 先 验 分 布 (Admissible-Range-Restricted 
Priors)， 也 可 以 较 大 程度 上 改善 DSEM 的 参数 估 
计 结 果 (McNeish, 2019)。 
4.1.2 ”模型 拟 合 及 比较 

DSEM 可 以 利用 偏差 信息 准则 (Deviance 
Information Criterion，DIC) 作 为 模型 拟 合 与 比较 
的 指标 (Asparouhov et al., 2018)。DIC 常用 于 贝 叶 
斯 估计 中 的 模型 比较 (Spiegelhalter et al., 2002), 
DIC 越 小 的 模型 拟 合 越 好 。 然 而 ， 当 DSEM 中 参 
数 数 量 较 多 、 模 型 过 于 复杂 时 ， 需 要 较 多 迭代 次 
数 才 能 得 到 较为 稳定 的 DIC (Asparouhov et al., 
2018)。 因此, Asparouhov A (2018) EX, 研究 者 
需要 通过 对 模型 设置 不 同 的 随机 种 子 并 多 次 运行 
模型 ， 以 确认 DIC 是 否 已 经 达到 稳定 的 取 值 。 
4.1.3 ”优势 及 不 足 

DSEM 在 密集 追踪 数据 的 分 析 上 存在 诸多 亮 
点 和 优势 ,首先 , DSEM 结合 多 种 模型 方法 的 优势 ， 
能 更 灵活 、 方 便 地 对 模型 中 多 个 变量 之 间 建 立 相 
关 关 系 (Hamaker et al., 2018; Zhou et al., 2020) 。 它 
不 仅 局 限于 对 两 变量 的 建 模 分 析 , 也 可 以 对 三 个 
及 以 上 变量 进行 建 模 。 不 管 是 显 变量 还 是 潜 变 量 、 
是 否 是 在 同一 时 间 节 点 测量 , DSEM 中 各 个 变量 
之 间 均 可 建立 回归 路 径 (Asparouhov et al., 2018)。 
其 次 , DSEM 能 方便 地 将 密集 追踪 数据 中 呈现 的 
特征 纳入 模型 (如 ， 考虑 数据 中 存在 的 动态 趋势 、 
周期 变化 等 ) 不 仅 可 以 考虑 变量 本 身 随时 间 的 变 


时 间 节 点 和 随机 效应 也 不 会 对 其 带 来 较 多 计算 
量 上 的 负担 (Asparouhov et al., 2018)。 前 人 模拟 研 
究 表明 ，DSEM 在 随机 效应 的 均值 估计 方面 表现 
十 分 优秀 (Schultzberg & Muthén, 2018)。 

然而 , DSEM 也 存在 一 些 不 足 与 局 限 ,在 对 于 
间隔 不 等 的 时 间 节 点 的 处 理 中 ,不 同时 间 间 隔 5 
的 选择 会 对 DSEM 水 平 1 参数 的 估计 结果 产生 一 
定 的 影响 (Asparouhov et al., 2018; Asparouhov & 
Muthen, 2019b)。Asparouhov 等 人 (2018) 在 DSEM 
的 基础 上 做 出 微调 ， 提 出 了 残 差 动 态 结构 方程 模 
型 (Residual Dynamic Structural Equation Model, 
RDSEM)， 仅 在 测量 模型 中 体现 自 回 归 和 交叉 沸 
后 的 效应 ,因此 其 结构 方程 模型 部 分 不 会 受到 不 
同时 间 间 隔 选 择 的 影响 (Asparouhov et al., 2018)。 
Asparouhov 和 Muthén (2019b) 的 模拟 研究 表明 ， 
在 观测 时 间 节 点 的 间隔 不 一 的 情况 下 , RDSEM 的 
模型 结果 相对 DSEM 来 说 更 加 稳健 ， 几 乎 不 会 受 
到 时 间 间 隔 的 影响 。 

总 的 来 说 DSEM 在 密集 追踪 数据 的 分 析 上 
存在 诸多 亮点 和 优势 (Hamaker et aL, 2018), 但 也 
仍 面临 许多 新 的 挑战 ， 明 待 更 多 心理 学 、 统 计 学 
等 领域 的 研究 者 共同 解决 。 

4.2” 自 下 而 上 的 方法 : 组 迭代 多 模型 估计 

本 质 上 来 说 , DSEM 和 传统 MLM 一 样 ,将 被 
试 视 为 整体 的 、 同 质 化 的 群体 ,认为 在 个 体 间 存 
在 相同 的 模型 路 径 , 不同 被 试 之 间 可 能 只 是 效应 
大 小 有 所 不 同 ( 这 些 随机 效应 可 以 利用 个 体 特异 
和 时 间 特 异 的 变量 进行 解释 )， 因 此 是 一 种 基于 整 
体 的 方法 (Piccirillo & Rodebaugh, 2019)。 然 而 ， 基 
于 整体 和 基于 个 人 的 分 析 方 法 可 能 得 到 截然 不 同 
的 结果 (Fisher et al., 2018; Turner & Hayes, 1996), 
特别 是 实际 中 个 体 特异 性 较 强 的 情况 下 ， 针 对 整 
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体 数 据 构 建 的 群体 模型 可 能 不 能 解释 群体 内 部 的 
任 一 个 体 (Miller & van Horn, 2007; Molenaar & 
Campbell, 2009)。 即使 患 有 同一 疾病 ,不 同 个 体 之 
间 的 具体 情况 、 变 化 模式 等 也 存在 较 大 的 个 体 差 
(Hayes et al., 2019)。 因 此 ， 近 年 来 越 来 越 多 学 
者 倡导 研究 者 开展 基于 个 体 的 研究 和 统计 分 析 
(Foster & Beltz, 2018; Hayes et al., 2019; Piccirillo 
& Rodebaugh, 2019; Wright & Woods, 2020)。 利 用 
基于 个 人 的 分 析 方 法 对 不 同时 间 、 不 同情 境 下 的 
动态 变化 过 程 进 行 建 模 和 分 析 , 可 以 得 到 更 精确 
的 、 更 具 针 对 性 的 结果 ， 从 而 增加 治疗 和 干预 对 
于 个 人 的 匹配 度 和 治疗 效用 (Hayes et al., 2019; 
Wright & Woods, 2020)。 

纯粹 基于 个 人 的 分 析 方 法 充分 考虑 了 个 体 间 
的 异 质 性 , 但 同时 也 受到 了 一 些 质 疑 。 如, (MSE 
于 个 人 建立 的 模型 可 能 不 够 可 靠 , 容易 受到 数据 
中 噪音 的 影响 ,从 而 得 到 过 拟 合 的 模型 (Gates & 
Molenaar，2012)， 缺 乏 一 定 的 普 适 性 (Spencer & 
Schöner, 2003)。 由 此 , Gates 和 Molenaar (2012) 在 
针对 个 人 的 模型 的 基础 上 ,结合 群体 共享 信息 的 
优势 ， 提出 GIMME 这 一 算法 。GIMME 最 初 是 作 
为 一 种 针对 功能 性 磁 共 振 成 像 (Functional Magnetic 
Resonance Imaging，fMRD 的 数据 分 析 方 法 提出 ， 
但 近年 来 也 被 较 多 地 应 用 于 行为 和 心理 构 念 方面 
密集 追踪 数据 的 分 析 (Foster & Beltz, 2018; Gates 
et al, 2020), 它 可 以 被 视 作 纯 粹 基于 个 体 的 方法 
和 基于 整体 的 方法 之 间 的 一 种 连接 (Beltz et al., 
2016; Wright & Woods, 2020)。 
424 方法 原理 

GIMME 是 一 种 动态 网 络 的 方法 ,每 个 变量 
均 被 视 为 网 络 中 的 一 个 节点 。 对 于 一 组 时 间 序 列 
观测 变量 (7) (p x 1), GIMME 的 公式 如 下 所 示 
(Beltz & Gates, 2017): 
nj (t) 7 CA; + AF) (t) +(G ; + DE. y, (t -1)+& (t) (12) 

其 中 ，4 (p xD) 表示 相同 时 间 点 各 个 变量 节 
点 之 间 的 边 (Contemporaneous Edges)， 可 以 理解 
为 变量 间 的 同时 效应 。 5B (p xp) 为 各 变量 节点 间 
一 阶 的 沾 后 效应 (First Order Lagged Edges)， 表 示 
在 时 间 点 t 时 ,变量 节点 与 上 一 时 间 点 的 测量 之 
间 的 关系 ， 可 能 同时 包括 自 回 归 和 交叉 滞后 的 效 
应 。&(1) (px1) 表 示 均 值 为 0 的 误差 项 。GIMME 
的 个 体 模型 是 基于 统一 结构 方程 模型 (Unified 
Structural Equation Modeling, USEM; Kim et al., 


2007) y. . fll DFM 不 同 , USEM 是 针对 显 变量 (而 
不 是 潜 变 量 ) 进 行 建 模 (Kim et al., 2007)。 在 USEM 
的 基础 上 ，GIMME 还 区 分 了 个 体 和 群体 的 效应 ， 
下 标 i 表 示 该 系数 为 第 i 个 被 试 独 有 ,而 上 标 g 表 
示 群 体 模型 中 的 效应 。 因 此 , 在 各 个 个 体 的 模型 
中 ,将 同时 包括 个 体 特 异 以 及 群体 共享 的 路 径 。 
简单 来 说 , GIMME 的 建 模 分 析 步 骤 大 致 包括 : (1) 
数据 预 处 理 ; (2) 基 于 大 多 数 个 体 模型 中 存在 的 路 
径 ， 自 下 而 上 地 建立 群体 的 模型 ; (3) 以 群体 模型 
为 初始 模型 建立 个 体 的 模型 (Beltz & Gates, 2017; 
Gates & Molenaar, 2012)。 

GIMME BUAN, RAW E (t) 与 时 间 相互 
独立 , 不 考虑 变量 随时 间 发 展 的 线性 趋势 (Beltz 
& Gates, 2017)。 因 此 , 人 研究 者 需要 将 原始 数据 对 
时 间 这 一 变量 进行 回归 分 析 ， 去 除 时 间 对 变量 的 
线性 影响 后 ， 再 将 得 到 的 残 差 项 纳入 GIMME 中 
进行 正式 的 分 析 (Beltz & Gates, 2017)。 另 外 ， 研 究 
者 也 可 以 将 时 间作 为 外 源 变 量 (Exogenous Variable) 
纳入 模型 从 而 去 除 各 变量 随时 间 变 化 的 趋势 
(Woods et al., 2020)。 

在 群体 模型 的 建立 中 , GIMME 将 初始 模型 作 
为 起 点 , 根据 修正 指数 (Modification Indices) 对 个 
体 模型 逐步 迭代 地 添加 回归 路 径 ， 只 有 显著 优化 
大 多 数 个 体 模型 拟 合 的 路 径 才 能 被 纳入 群体 模型 
中 (Gates & Molenaar, 2012)。 接 着 , GIMME 以 群 
体 模 型 作为 迭代 的 初始 模型 ,对 各 个 个 体 模 型 逐 
步 添加 显著 优化 其 模型 拟 合 的 路 径 ， 得 到 各 个 个 
体 的 模型 。 关 于 GIMME 建立 群体 和 个 体 模型 的 
具体 算法 可 见 Gates 和 Molenaar (2012) 的 文章 。 

一 般 来 说 ，GIMME 的 一 阶 延 迟 模型 (包括 同 
时 效应 和 谐 后 一 个 时 间 点 的 效应 ) 已 经 能 较 好 地 
解释 密集 追踪 数据 中 的 时 间 依 赖 性 ,但 还 是 有 必 
要 对 每 个 个 体 的 模型 进行 后 验 模 型 验证 (A 
Posteriori Model Validatiom)， 以 查看 是 和 否 需要 考 
虑 更 高 阶 的 模型 (Beltz & Gates, 2017; Beltz & 
Molenaar, 2015)。 这 种 后 验 模型 的 验证 可 以 通过 检 
查 残 差 项 EO) 是 否 为 白 噪 音 来 实现 (Beltz & 
Gates, 2017). Beltz 和 Molenaar (2015) 详 细 介绍 J 
GIMME 中 后 验 模型 验证 的 方法 和 原理 ， 并 提供 
了 相应 的 LISREL 代码 。 

在 GIMME 的 基础 上 ， 后 续 较 多 研究 也 提出 
了 其 拓展 模型 ， 如 Gates 等 人 (2017) 提 出 的 考虑 子 
群 组 的 GIMME (Subgrouping Group Iterative Multiple 
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Model Estimation，S-GIMME)，Beltz 和 Molenaar 
(2016) 提 出 的 多 解决 方案 的 GIMME (GIMME for 
Multiple Solutions, GIMME-MS), 以 及 Gates 等 人 
(2020) 提 出 的 基于 潜 变 量 进行 建 模 的 GIMME 
(Latent Variable GIMME, LV-GIMME). 在 GIMME 
的 基础 上 ，LV-GIMME 利用 带 有 工具 变量 
(Instrumental Variables) 的 两 阶段 最 小 二 乘法 估计 
(Two-Stage Least Squares Estimation, MIIV-2SLS), 
同时 允许 测量 模型 和 结构 模型 部 分 在 不 同 个 体 间 


在 诸多 优势 , 但 它 同 样 存在 一 些 局 限 。 首 先 
GIMME 本 质 上 是 一 种 数据 驱动 的 建 模 方法 ， 可 
能 出 现 模 型 过 拟 合 (Overfitting) 等 问题 (Beltz & 
Gates, 2017)。 其 次 , GIMME 主要 考虑 变量 间 共 变 
(Covary) 的 关系 及 机 制 ， 无 法 将 几乎 不 随时 间 变 
化 的 变量 纳入 模型 进行 分 析 (Beltz & Gates, 
2017)。 为 了 探究 个 体 水 平 上 不 随时 间 变 化 的 协 变 
量 的 作用 ,人 研究 者 需要 额外 对 个 体 模型 中 得 到 的 
参数 估计 进行 进一步 分 析 。 如 , Kim 等 人 (2007) 利 


存在 差异 。 模 拟 研究 的 结果 表明 ,， 相 比 于 对 条 目 
总 分 直接 进行 建 模 ,，LV-GIMME 能 更 有 效 地 检测 
出 模型 中 存在 的 效应 (Gates et al., 2020)。 
4.2.2 ”优势 及 不 足 

总 的 来 说 , GIMME 在 密集 追踪 数据 的 分 析 上 
具有 其 独特 的 优势 。 与 基于 整体 的 分 析 方 法 相 比 ， 
GIMME 的 优势 主要 体现 在 : (1) 对 于 群体 模型 的 
建立 是 一 个 自 下 而 上 的 过 程 ， 其 路 径 的 纳入 与 否 
取决 于 该 路 径 能 否 显著 改善 一 定 比例 的 单个 个 体 
模型 (Gates & Molenaar，2012)， 可 以 尽 可 能 减少 
个 体 特异 信息 的 损失 ,更 适用 于 探究 变量 间 共 变 
模式 的 个 体 差 异 (Beltz et al., 2016); (2) 人 允许 各 个 
个 体 的 模型 结构 之 间 存 在 异 质 性 ， 并且 能 提供 任 
一 个 体 相对 于 整个 群体 的 差异 方面 的 信息 ， 从 而 
更 好 地 理解 个 体 模型 的 特征 和 模式 ， 这 对 于 临床 
上 针对 个 人 的 精准 治疗 来 说 具有 重要 的 意义 
(Wright & Woods, 2020); (3) 对 于 被 试 数量 的 要 求 
BUD, 在 只 有 10 名 被 试 的 情况 下 就 已 经 可 以 较 好 
地 对 群体 模型 的 路 径 进 行 估计 (Gates & Molenaar, 
2012; Henry et al., 2019)。 而 与 纯粹 基于 个 体 的 方 
法 相 比 , GIMME 也 具有 其 独特 的 优势 : 它 在 考虑 
个 体 特异 性 的 同时 大 大 提高 了 结果 的 普 适 性 , 一 
定 程 度 上 降低 了 模型 过 拟 合 的 风险 (Beltz & Gates, 
2017)。 另 外 , GIMME 在 软件 操作 上 较为 简便 ， 且 
能 自动 输出 每 个 个 体 模型 以 及 群体 模型 的 各 项 拟 
合 指标 、 路 径 系数 等 ， 为 研究 者 基于 个 人 的 精细 
化 研究 提供 了 丰富 的 信息 。GIMME 还 可 以 较为 直 
观 地 对 结果 进行 呈现 ,更 易于 理解 与 解释 ,适用 
于 动态 系统 中 存在 较 多 相互 影响 的 变量 的 情况 
(Beltz et al., 2016), GIMME 不 仅 可 以 帮助 研究 者 
从 宏观 的 视角 理解 变量 间 的 关系 ,还 可 以 输出 每 
个 个 体 的 模型 参数 估计 结果 ,大 大 增加 治疗 效 月 
和 干预 指导 的 精确 性 。 
虽然 GIMME 在 对 密集 追踪 数据 的 分 析 上 存 


ag 


用 一 般 线 性 模型 (General Linear Model), 探究 个 
体 层 面 的 变量 (性 别 、 年 龄 、 智 商 及 教育 ) 对 USEM 
中 得 到 的 个 体 模型 路 径 的 影响 。 又 如 ，Beltz 等 人 
(2013) 利 用 7 检验 ,比较 USEM 得 到 的 个 体 模 型 路 
径 中 , 各 变量 间 达 到 显著 的 回归 路 径 的 数量 是 否 
存在 显著 的 性 别 差异 。 另 外 , GIMME 假设 认为 变 
量 间 的 效应 相对 时 间 静 止 ， 不 考虑 数据 中 可 能 存 
在 的 时 变性 质 (Time-Varying Nature; Beltz & Gates, 
2017; Foster & Beltz, 2018), 但 实际 研究 中 可 能 出 
现 变量 间 关 系 随时 间 动 态 变 化 的 情况 。 最 后 ， 
GIMME 要 求 各 个 测量 时 间 节 点 之 间 必 须 间 隔 相 
等 ,但 在 许多 密集 追踪 的 研究 中 ， 观 测 节点 的 时 
间 间 隔 常常 是 不 等 的 (Foster & Beltz, 2018; Kleiman 
et al., 2019)。 另外 ,虽然 GIMME 对 于 个 体 数 量 上 
几乎 没有 限制 , 但 它 要 求 对 于 每 个 个 体 进行 较 多 
时 间 节 点 的 测量 。 
43 ”实证 演示 及 方法 比较 

DSEM 已 经 可 以 在 Mplus 上 实现 ， 且 代码 语 
句 较 为 简单 ， 这 为 研究 者 们 利用 DSEM 进行 密集 
追踪 数据 的 动态 分 析 提 供 了 很 大 的 便利 (McNeish 
& Hamaker, 2019)， 近 年 来 涌现 出 较 多 关于 DSEM 
的 应 用 研究 (Armstrong et al., 2019; Hamaker et al., 
2018; Racine et al., 2019; Zhang & Zheng, 2019)。 
而 GIMME 已 有 与 其 同名 的 R 包 (gimme; Lane 
etal, 2020), 也 可 以 利用 LISREL 以 及 MATLAB 
中 的 工具 箱 进 行 建 模 ， 近 年 来 在 应 用 领域 也 出 现 
较 多 研究 利用 GIMME (Beltz & Gates, 2017; 
Bouwmans et al.，2018) 或 GIMME 的 扩展 模型 进 
行 密集 追踪 数据 的 分 析 ， 如 考虑 子 群 体 模式 的 
S-GIMME (Lane et al., 2019)， 以 及 可 以 同时 输出 
多 个 群体 模型 方案 的 GIMME-MS (Beltz et al., 
2016)。 

为 了 进一步 加 深 读者 对 于 DSEM 和 GIMME 
实际 建 模 过 程 的 理解 、 突 出 二 者 相 比 传统 MLM 
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的 建 模 优 势 ， 本 节 将 利用 Mplus 8.4 和 R 包 gimme， 
对 McNeish 和 Hamaker (2019) 提 供 的 公开 数据 分 
别 进行 MLM, DSEM 和 GIMME 的 建 模 分 析 ( 具 
体 代 码 见 网 络 版 附录 )。 简 单 来 说 ，McNeish 和 
Hamaker (2019) 对 100 名 被 试 在 50 个 时 间 点 的 抑 
郁 水 平 (Depression) 和 吸烟 欲望 (Urge to Smoke) 进 
行 密集 追踪 , 利用 DSEM 一 系列 不 同形 式 的 模型 
探讨 抑郁 水 平和 吸烟 欲望 之 间 的 动态 关系 ， 以 及 
个 体 特异 的 变量 (工作 压力 和 家 庭 压力 ) 在 其 中 所 
起 的 作用 。 
4.3.1 ”传统 多 层 线性 模型 的 实证 演示 

McNeish 和 Hamaker (2019) 的 数据 中 同时 包 
括 个 体 层 面 和 时 间 点 层面 的 变量 测量 ， 且 因 变 量 
吸烟 欲望 的 组 内 相关 系数 (Intraclass Correlation, 
ICC) 较 大 (CC = 0.387)， 因 此 可 以 利用 MLM 进行 
数据 分 析 。 假 设 吸烟 欲望 的 基线 水 平 、 随 时 间 变 
化 的 趋势 以 及 抑郁 对 吸烟 欲望 的 影响 在 不 同 个 体 
间 存 在 较 大 的 差异 ， 因 此 本 例 将 建立 带 随 机 截 距 
和 随机 斜率 的 MLM 模型 (对 应 Mplus 代码 见 网 络 
版 附录 A)。 为 了 更 好 地 解释 结果 ,我们 对 水 平 1 
的 预测 变量 (抑郁 水 平 ) 进 行 个 体 均 值 中 心 化 ,根据 
3.1 节 中 对 于 MLM 方 法 原理 的 介绍 ,在 水 平 1, 个 
体 i 在 测量 时 间 点 t 的 吸烟 欲望 不 可 以 表示 为 : 

Yi, = Boi + PiX n — X;) + Boi; + Ei (13) 

Hp, XQ X, 分别 表示 个 体 i 在 第 ieee! 
间 点 的 抑郁 水 平 以 及 个 体 i 的 平均 抑郁 水 平 ， 
表示 测量 时 间 。 我 们 假设 水 平 1 CINE RE 
在 不 同 个 体 间 存 在 变异 ， 且 这 一 变异 可 能 可 以 被 
个 体 本 身 来 自 工 作 和 家 庭 的 压力 解释 。 因 此 , 在 
水 平 2, MLM 的 模型 表达 式 如 下 : 


| 


T 


Boi = Yoo + 10W i, i Yao zit, i + Uoi (14) 
Bi = Yoi +71 Wate, i + Yo gee, ity (15) 
Bai = Yor + YoW ye i + VW gene, o (16) 

Kp, Wip Wag i 均 为 水 平 2 的 预测 变量 ,分 
别 表示 个 体 的 工作 压力 和 家 庭 压力 。 MLM 对 应 的 
模型 图 如 图 1 所 示 。 

Mplus 8.4 默认 对 MLM 采用 稳健 极 大 似 然 估 
计 (Robust Maximum Likelihood Estimator, MLR), 
因此 仅 允 许 输出 非 标准 化 的 结果 (Muthén & 
Muthén, 1998-2017). 如 果 需 要 输出 标准 化 的 结 
可 以 采用 贝 叶 斯 估计 的 方法 。 结 果 显 示 ， THE 
力 和 家 庭 压力 越 高 的 个 体 往往 有 更 强烈 的 吸烟 欲 
FA( 4) = 0.458, 79 = 0.286, ps < 0.001)。 个 体 的 吸 
烟 欲 望 随 时 间 的 推移 而 逐渐 增加 (yo = 0.157, p < 
0.001), 但 这 种 随时 间 变 化 的 趋势 在 个 体 间 的 随 
机 效应 不 能 被 个 体 本 身 的 工作 压力 或 家 庭 压力 解 
释 (ps > 0.05)。 抑 郁 水 平 越 高 时 , 个 体 的 吸烟 欲望 
也 越 强 烈 (y01= 0.741, p < 0.001), 且 个 体 的 工作 
压力 和 家 庭 压力 越 大 ,该 个 体 抑郁 水 平 对 吸烟 欲 
望 的 影响 也 越 大 (%1= 0.270, p = 0.001; y= 
0.334, p < 0.001)。 总 之 , 个 体 本 身 的 工作 压力 和 
家 庭 压力 会 提高 个 体 吸烟 欲望 的 基线 水 平 ， 且 增 
加 其 抑郁 水 平 对 吸烟 欲望 的 影响 。 
43.2 ”动态 结构 方程 模型 的 实证 演示 

与 MLM 相 比 , DSEM 可 以 方便 地 对 密集 追踪 
数据 中 的 动态 效应 进行 建 模 分 析 。 本 例 仅 考虑 了 
个 体 特 异 的 预测 变量 ， 因 此 采用 两 层 的 DSEM. 
McNeish 和 Hamaker (2019) 的 数据 中 ， 吸 烟 和 欲望 
整体 呈现 出 随时 间 线 性 增长 的 趋势 ， 因 此 需要 对 
数据 进行 去 时 间 趋 势 (Detrending)。Asparouhov 
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吸烟 欲望， z 


对 吸烟 欲望 的 影响 


TE: 圆圈 内 变量 为 潜 变量 , 方 框 内 变量 为 显 变 下 标 i 表示 第 i 个 个 体 的 测量 , 1 表示 第 1 个 时 间 点 的 观测 数据 。 
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(2018) 的 模拟 研究 表明 , DSEM 在 同时 进行 去 时 间 
趋势 和 自 回归 效应 的 分 析 时 会 对 方差 的 估计 产生 
有 差 ， 因 此 建议 采用 DSEM 的 变 式 一 -RDSEM i 
行 分 析 。RDSEM 与 DSEM 的 区 别 主要 在 于 
RDSEM 分 别 对 去 时 间 趋 势 与 自 回 归 和 交叉 滞后 
效应 进行 建 模 ， 即 仅 在 测量 模型 部 分 而 不 是 结 
构 方 程 模型 的 部 分 考虑 自 回 归 和 交叉 滞后 效应 。 
RDSEM 是 对 DSEM 微调 后 的 提出 的 一 种 变 式 
(Asparouhov et al., 2018)， 因 此 对 RDSEM 的 介 
绍 同样 可 以 较 好 地 展示 DSEM 的 分 析 过 程 及 
思路 。 

根据 理论 假设 ,在 水 平 1, 吸烟 欲望 和 抑郁 
水 平 存在 自 回归 ( ga gag Dag agi ) ABE OCT Jr 
的 效应 ， 即 前 一 个 时 间 点 的 抑郁 水 平 会 对 后 一 时 
PEJ BD IR ARI RK EE PA HE MA (Dyson anya )， 而 前 一 时 
间 点 的 吸烟 欲望 也 可 能 会 对 后 一 时 间 点 的 抑郁 水 
FETE Dig ay oui )。 对 应 模型 公式 为 : 


(w) 

bon £ Êy- i T, + 5 吸烟 , it (18) 
Yir, it Pini; 6 抑郁 ,it 
bail Dr Ra 


Ot qug; Pan- i | eiim i, 1-1 


hn "| (19) 


6 抑郁 , it 


Ogg, i 

与 4.1.1 WHR DSEM 的 公式 相 比 , RDSEM 
仅 在 测量 模型 部 分 考虑 残 差 项 5 的 自 回 归 和 交 
SC tit Ja AY BOM oS gus; = Var (guam i) Tig, i = 
var Sip, )， 则 在 水 平 2， 个 体 的 吸烟 欲望 和 抑郁 
的 平均 水 平 ( Hgm o Lg; )、 各 个 回归 系数 
C Pig pa is Dogg, i Pv io Dm, i Pc i+ 
Brn- ) A BAe HIT 25 Cru js Tans, ) 均 可 
能 受到 个 体 的 工作 压力 和 家 庭 压力 的 影响 。 在 考 
虑 残 差 项 方差 的 随机 效应 时 ， 为 了 保证 残 差 项 方 
差 为 正 数 ，DSEM 将 取 对 数 之 后 的 残 差 项 方差 


| | xo M 
ead "um. i H (17) 75g, ; 和 Tg 纳入 水 平 2 的 模型 。 因 此 , 水 平 2 
Yes] Manin | Yima 模型 的 公式 为 : 
Hay, i Fu = /工作 ,吸烟 7 家 庭 ,吸烟 吸烟 ,i 
吸烟 
NI 7 工作 , 抑郁 /家 庭 , 抑郁 2 抑郁, 
抑郁 Vig i 
Day gup ; Tipi /工作 ,吸烟 -吸烟 /家 庭 ,吸烟 -吸烟 | | 2 吸烟 -吸烟 
Pomi | | Vn 工作 ,吸烟 -抑郁 7 家 庭 ,吸烟 -抑郁 | | Uomo i 
O ua num ; 7 抑郁 -吸烟 /工作 , 抑郁 -吸烟 /家 庭 , 抑郁 -吸烟 抑郁 -吸烟 ,i 
P MM | Wt, tug; * Q0) 
Os ang i 抑郁 -抑郁 7 工作 , 抑郁 -抑郁 7 家 庭 , 抑郁 - 抑 邦 | | 抑郁 -抑郁 ,i 
Bogman | | 758-8 工作, 吸烟- 时间 7 家庭 ,吸烟 -时 间 | | "gn al i 
7 抑郁 -时 间 
Py utin; 7 工作 , 抑郁- 时间 7 家庭, 抑 邦 -时 间 | | Mgmt, i 
1 Nios ny) iü 
Ogg i ;) T Ef, log ayy) T x, log zig) log zu, ;) 
log(z ) | os) | M 
= 8 LM L Pr HF, log (yr ) | | is, logzgg) | | lg) | 


其 中 ,公式 (17) 对 应 表示 变量 变异 来 源 的 分 
解 ， 公式 (18)~(19) 对 应 水 平 1 的 模型 ， 公 式 (20) 对 
应 水 平 2 的 模型 (对 应 Mplus 代码 见 网 络 版 附录 
A)。McNeish 和 Hamaker (2019) 的 数据 仅 提供 了 
个 体 特异 预测 变量 的 数据 ,因此 本 例 仅 假设 水 平 
1 的 各 项 系数 在 水 平 2 受到 个 体 特异 预测 变量 的 
影响 ; 但 如 果 数 据 信 息 充 足 ， 我 们 也 可 以 类 似 地 
在 水 平 3 考虑 时 间 特 异 的 预测 变量 (如 ,特定 的 工 
作 日 ) 产 生 的 影响 ,如 可 以 将 残 差 项 方差 za ; 的 
随机 效应 抽象 地 表示 为 logg, i) = 52,4 + 83,45 DA 
而 将 其 分 解 为 个 体 和 时 间 特 异 的 随机 效应 。 在 本 


例 中 ,对 于 第 i 个 个 体 ， 其 模型 图 如 图 2 所 示 。 
DSEM 和 RDSEM 均 可 以 同时 输出 标准 化 和 
非 标准 化 的 结果 ,， 且 标准 化 结果 还 在 水 平 1 输出 
了 不 同 个 体 中 各 项 标准 化 系数 的 平均 结果 (Within- 
Level Standardized Estimates Averaged Over Clusters)。 
为 了 更 好 地 对 吸烟 欲望 和 抑郁 水 平 之 间 的 交叉 滞 
后 效应 进行 解读 , 我们 采用 的 是 Mplus 输出 的 标 
准 化 结果 ， 其 标准 化 的 具体 过 程 及 原理 见 
Schuurman 等 人 (2016) 的 文章 。 水平 1 标准 化 系数 
在 不 同 个 体 中 的 平均 结果 显示 ,吸烟 欲望 的 自 回 
归 效 应 达到 显著 ( By guy = 0.177, CI = [0.142， 
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变量 分 解 


HRM, i 


图 2 ”吸烟 欲望 和 抑郁 水 平 的 动态 关系 


ik. 圆圈 内 变量 为 潜 变量 ， 方 框 内 变量 为 显 变量 。(W) 表 示 个 体内 的 估计 , “吸烟 " 和 “抑郁 * 分 别 表示 吸烟 欲望 和 抑郁 水 平 。 
Dai 和 更 :分别 表示 吸烟 欲望 和 抑郁 水 平 的 自 回归 效应 ，Bx; 和 Da; 分 别 表示 上 一 时 间 点 的 吸烟 欲望 对 当前 抑郁 水 平 上 


一 时 间 点 的 抑郁 水 平 对 当前 吸烟 欲望 的 交叉 沾 后 效应 。 p 表示 均值 ，6w ，， 和 Siggy,» 分 别 表示 吸烟 欲望 和 抑郁 水 平 的 残 差 


The mi Mm, NAO, MS, ,的 方差 。 为 了 更 清晰 地 呈现 模型 假设 , 模型 图 中 未 呈现 与 去 时 间 趋 势 有 关 的 路 径 。 


0.211])， 前 一 时 间 点 的 吸烟 欲望 可 以 显著 预测 后 
一 时 间 点 的 吸烟 欲望 。 另 外 ,被 试 的 吸烟 欲望 整 
体 随 时 间 的 推移 而 逐渐 增加 ( Byyg_wyy = 0-721, 
CI = [0.704, 0.738])。 前 一 时 间 点 的 抑郁 水 平 对 后 
一 时 间 点 吸烟 欲望 的 影响 没有 达到 显著 , 但 前 一 
时 间 点 的 吸烟 欲望 可 以 降低 后 一 时 间 点 的 抑郁 水 
(Du ug = —0.069, CI = [-0.102, -0.038])。 在 
水 平 2, 工作 压力 和 家 庭 压 力 越 大 的 个 体 ， 其 吸 
AKERI tema = 0.362, CI = [0.190, 0.515]; 
Yagrgug = 0.220, CI = [0.055, 0.371])， 吸烟 欲望 
的 自 回归 效应 越 强 (7 人 ag ug = 0.423, CI = 
[0.220, 0.594]; vag 吸烟 -吸烟 = 0.212, CI = [0.016, 
0.388]), 吸烟 欲望 的 残 差 项 方差 越 大 ( yw, tong) 
= 0.355, CI = [0.227, 0.474]; yi, ogra) = 0-211, 
CI = [0.076, 0.334])， 因 此 其 吸烟 欲望 的 残 差 的 动 
态 波 动 幅度 也 越 大。RDSEM 的 结果 进一步 表明 ， 
个 体 的 吸烟 欲望 存在 较 强 的 惯性 ( 自 回归 效应 )， 
且 个 体 的 工作 压力 和 家 庭 压力 越 大 ,这 种 惯性 也 


[S 


点 的 吸烟 欲望 可 以 降低 后 一 时 间 点 的 抑郁 水 平 。 
现实 生活 中 ,变量 间 的 相互 作用 很 可 能 存在 复杂 
的 双向 反馈 ， 而 DSEM 允许 在 同一 模型 中 建立 双 
向 的 作用 路 径 ， 这 可 能 可 以 为 研究 者 带 来 更 多 
新 的 研究 思路 和 启发 。 
4.3.3 ”组 迭代 多 模型 估计 的 实证 演示 

本 节 将 利用 gimme 这 一 R 包 进 行 GIMME 的 
建 模 分 析 。 在 对 数据 进行 分 析 之 前 ,首先 需要 去 
除数 据 中 随时 间 线 性 变化 的 趋势 。 研 究 者 可 以 将 
原始 数据 对 时 间 进 行 回归 后 得 到 的 残 差 项 纳 人 
GIMME 中 进行 正式 的 分 析 ， 也 可 以 直接 将 时 间 
作为 外 源 变 量 纳入 模型 ， 从 而 去 除 各 变量 随时 间 
变化 的 趋势 (Woods et al., 2020)。 为 了 更 直观 地 是 
现 去 时 间 趋 势 的 效果 ,本 例 采用 前 一 种 去 时 间 趋 
势 的 方法 。 以 第 一 个 被 试 的 吸烟 欲望 为 例 , 图 3a 
为 吸烟 欲望 的 原始 测量 (URGE) 随 时 间 变 化 的 总 
体 趋势 ， 随 着 时 间 的 推移 ,吸烟 欲望 越 强烈 。 通 过 
建立 时 间 对 吸烟 欲望 的 回归 模型 得 到 的 残 差 来 去 


越 强 。 我 们 在 进行 MLM 的 分 析 中 发 现 抑郁 水 平 
对 吸烟 欲望 存在 正 向 预测 的 作用 , 而 RDSEM 的 
结果 体现 了 相反 方向 的 另 一 种 可 能 ， 即 前 一 时 间 


除数 据 中 随时 间 线 性 变化 的 趋势 后 ， 其 结果 如 图 
3b 所 示 。 可 以 看 出 ， 此 时 吸烟 欲望 在 不 同时 间 点 
的 测量 呈 平 稳 的 趋势 。 
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图 3 吸烟 欲望 随时 间 发 生 的 变化 


去 除 各 变量 受到 的 来 自 时 间 的 影响 后 ， 将 其 
纳入 GIMME 进行 分 析 。 对 于 每 个 被 试 i， 其 吸烟 
欲望 和 抑郁 水 平 的 时 间 序 列 观测 变量 可 以 表示 为 : 


Naep, i(t) Tacp, it) 
Murge i(t -D Surge, it) 
D + @%. : 21 
(@ + zar ard (21) 


HP, Mugei lt) 和 Macy (t) 分 别 表示 被 试 i 在 
时 间 点 + 的 吸烟 欲望 及 抑郁 水 平 。4 (2 x 2) 表 示 吸 
烟 欲 望 和 抑郁 水 平 的 同时 效应 矩阵 ，@G(2 x 2) 为 吸 
烟 和 欲望 与 抑郁 水 平 之 间 的 一 阶 滞后 效应 (First 
Order Lagged Edges)， 包 括 自 回归 和 交叉 滞后 效 
应 ,下 标 i 表示 该 系数 为 第 i 个 被 试 独 有 ， 而 上 标 
g 表示 为 群体 模型 中 的 效应 ,在 利用 GIMME 进行 
数据 分 析 时 ,研究 者 可 以 利用 已 有 的 经 验 对 群体 
模型 中 可 能 存在 的 路 径 做 出 假设 ， 并 将 其 作为 初 
始 迭 代 的 模型 (如 ,在 gimmeSEM 图 数 中 声明 
paths = urge~deplag， 则 表示 初始 模型 中 存在 前 
一 时 间 点 的 抑郁 水 平 指向 当前 吸烟 欲望 的 回归 路 
径 )。 知 理论 依据 不 足 ， 也 可 以 直接 从 空 模型 出 发 ， 
以 数据 为 驱动 进行 探索 式 的 建 模 (Beltz & Gates, 
2017; Gates & Molenaar, 2012)。 本 文 将 利用 数据 
驱动 、 探 索 式 的 方式 进行 演示 。 

在 准备 好 数据 并 运行 函数 后 ，gimme 将 初始 
模型 作为 起 点 ， 自 下 而 上 地 建立 群体 模型 。 具 体 
KUL, gimme 首先 根据 修正 指数 对 每 个 个 体 模型 
逐步 迭代 地 添加 回归 路 径 ， 再 以 个 体 模型 中 的 路 
径 为 基础 建立 群体 模型 。 为 了 找到 在 大 多 数 个 体 
模型 中 较为 稳定 存在 的 路 径 和 效应 ,该 路 径 只 有 
能 显著 优化 大 多 数 个 体 模型 的 拟 合 时 (大 多 数 个 
体 模型 中 ,该 路 径 对 应 的 MI 在 0.01 的 水 平 上 达 
到 显著 ), 才能 被 纳入 群体 模型 中 (Gates & Molenaar, 


iE: 图 3a 为 吸烟 欲望 的 原始 测量 随时 间 变 化 的 总 体 趋势 ; 图 3b 表示 去 除数 据 中 随时 间 线 性 变化 的 趋势 之 后 的 结果 。 


2012)。 本 文 按 照 默 认 的 设置 , 将 75% 设 为 可 以 被 
纳入 群体 模型 的 临界 值 ， 研 究 者 也 可 以 根据 具体 
需要 进行 修改 。 基 于 群体 模型 , gimme 再 对 个 体 模 
型 逐个 进行 迭代 路 径 搜 索 。 在 完成 所 有 模型 的 路 径 
搜索 后 ，gimme 将 输出 各 个 模型 的 路 径 图 及 详细 
的 统计 结果 。 图 4 即 gimme 输出 的 群体 模型 图 。 和 群 
体 模型 结果 显示 , 超过 75% 的 个 体 模型 中 , 吸烟 欲 
望 和 抑郁 水 平均 存在 自 回 归 的 效应 (purge yurgelag = 
0.196, Baep>deplag = 一 0.049)。 而 其 他 路 径 均 没 有 在 
超过 75% 的 个 体 模型 中 出 现 ， 这 一 结果 也 体现 了 
变量 间 关 系 的 个 体 差异 。 


> 


图 4 群体 模型 
注 : urge: 吸烟 欲望 ; dep: 抑郁 水 平 。 黑 线 表示 群体 模型 中 
的 路 径 ， 灰 线 表示 个 体 模 型 中 存在 的 路 径 ， 实 线 代表 同时 


效应 ,虚线 表示 滞后 效应 (包括 自 回 归 和 交叉 澡 后 效应 )。 线 
条 粗细 表示 模型 中 存在 该 路 径 的 个 体 数 量 的 多 少 。 


图 5 为 随机 选取 的 个 体 模型 示例 。 总 体 来 看 ， 
不 同 个 体 在 吸烟 欲望 和 抑郁 水 平 之 间 的 动态 关系 
存在 较 大 的 个 体 差 异 。 如 ， 对 于 个 体 模 型 中 吸烟 
欲望 的 自 回归 效应 Purge ugelag， 较 多 个 体 在 前 一 
个 时 间 点 较 高 的 吸烟 欲望 均 会 预测 下 一 个 时 间 点 
更 高 的 吸烟 欲望 (47%)， 但 也 有 较 多 个 体 模型 中 
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图 5 个 体 模 型 (随机 示例 ) 


注 :urge: 吸烟 欲望 ; dep: 抑郁 水 平 。 红 线 表 示 正 向 的 效应 ， 蓝 线 表 示 负 向 的 效应 ， 实 线 代 表 同 时 效应 ， 虚 线 表 示 浪 后 效应 ( 包 


括 自 回 归 和 交叉 滞后 效应 )。 


二 者 不 存在 预测 作用 (47%)， 甚 至 可 能 存在 负 向 
的 自 回 归 路 径 (6%)。 另 外 , 在 9 个 个 体 模 型 中 ， 吸 
烟 欲 望 可 以 预测 下 一 时 间 点 的 抑郁 水 平 ， 其 中 4 
人 为 负 向 预测 , 5 人 为 正 向 预测 。 
4.3.4 “小结 

本 节 分 别 利用 MLM、RDSEM 和 GIMME if 
行 实证 演示 和 分 析 。RDSEM 以 DSEM 作为 基础 ， 
因此 本 节 也 较 好 展示 了 DSEM 的 分 析 过 程 及 思 
路 -MLM 和 DSEM 均 为 自 上 而 下 的 建 模 方 法 ， 相 
比 MLM, DSEM 在 分 析 密 集 追 踪 数 据 时 存在 诸多 
优势 : (1) 可 以 方便 地 对 变量 间 的 自 回 归 和 交叉 滞 
后 效应 进行 分 析 ( 包 括 相反 方向 的 交叉 滞后 效应 ); 
(2) 可 以 同时 考虑 多 个 因 变量 ; (3) 对 模型 随机 效应 
的 数量 几乎 没有 限制 ; (4) 可 以 在 水 平 1 输出 不 同 
个 体 中 各 项 标准 化 系数 的 平均 结果 ; (5) 采 用 潜 变 
量 中 心 化 ; (9) 可 以 考虑 残 差 项 方差 的 随机 效应 。 
如 果 研 究 者 希望 更 加 系统 、 深 入 地 学 习 DSEM 在 
Mplus 中 的 代码 实现 , 可 以 参考 McNeish 和 Hamaker 
(2019) 的 文章 。 

与 MLM 和 DSEM 不 同 , GIMME 是 一 种 典型 
的 自 下 而 上 的 建 模 方法 。 与 DSEM 相 比 , GIMME 
的 优势 主要 体现 在 它 可 以 单独 输出 各 个 个 体 模型 
的 结果 ， 人 允许 不 同 个 体 模型 的 结构 存在 异 质 性 ， 


从 而 更 有 针对 性 地 呈现 每 个 个 体 中 变量 相互 作用 
的 模式 。 然 而 , GIMME 不 能 将 不 随时 间 变 化 的 变 
量 纳入 模型 ， 因 此 不 能 方便 地 考虑 个 体 特异 的 变 
量 (如 工作 压力 、 家 庭 压力 ) 对 个 体内 部 变量 间作 用 
关系 的 影响 。 

44 其 他 模型 

除了 以 上 介绍 的 模型 方法 , 在 密集 追踪 数据 
中 也 存在 其 他 模型 方法 的 应 用 。 如 ,考虑 到 心理 
学 研究 中 常常 出 现 潜 变 量 的 情况 , 研究 者 们 将 
MLM 和 SEM 相 结 合 ,提出 多 层 结构 方程 模型 
(Multilevel Structural Equational Model, MSEM; 
Gottfredson et al., 2009; Johnson et al., 2005), 与 
MLM 类 似 , MSEM 在 探究 个 体内 部 的 动态 过 程 以 
及 个 体 特 异性 等 方面 也 存在 较 大 的 局 限 性 。 如 研 
究 者 希望 更 深入 地 了 解 MSEM， 可 以 参阅 Rabe- 
Hesketh 等 (2007), Preacher 等 (2010) 以 及 Holtmann 
等 (2016) 的 文章 。 

多 层 VAR 模型 (Multilevel Vector Autoregressive 
Model, mIVAR; Bringmann et al., 2015) 则 是 在 
VAR 模型 的 基础 上 ,进一步 允许 对 多 个 个 体 的 密 
集 追 踪 数 据 进 行 分 析 ， 也 是 目前 较为 常见 的 用 于 
密集 追踪 的 方法 之 一 (Bringmann et al., 2016; 
Jongerling et al., 2015; Piccirillo & Rodebaugh, 
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2019)。 特 别 是 其 与 网 络 分 析 相 结 合 ， 以 及 其 R 包 
mlVAR 的 提出 (Bringmann et al., 2013; Epskamp et 
aL, 2016), 使 得 mIVAR 的 应 用 相对 更 加 简便 ， 且 
同样 可 以 输出 与 GIMME 类 似 的 模型 网 络 ， 便 于 
在 临床 上 识别 症状 网 络 中 的 症状 中 心 (Symptom 
Centrality; Bringmann et al.，2015) 。 然 而 ， 与 
GIMME 不 同 , mIVAR 本 质 上 采用 的 是 基于 整体 的 
建 模 思想 。mlVAR 具有 测量 时 间 点 的 数量 要 求 较 
低 、 适 用 于 高 维 数据 分 析 的 优点 , 但 它 只 能 进行 
显 变 量 的 建 模 ， 且 难以 对 复杂 模型 进行 估计 
(Bringmann et al., 2015; Epskamp et al., 2018)。 由 
于 篇 幅 有 限 ， 本 文 仅 对 mlVAR 进行 简单 介绍 ,更 
详细 的 介绍 可 参见 Bringmann 等 人 (2015) 的 文章 。 

另外 ,考虑 到 随 着 时 间 的 推移 ， 变 量 动态 变 
化 过 程 中 可 能 存在 一 定 的 时 变 效应 (Time-Varying 
Effect), Bringmann 等 人 (2018) 在 VAR 的 基础 上 进 
一 步 将 时 变 效应 纳入 考虑 ， 提 出 时 变 向 量 自 回归 
模型 (Time-Varying Vector Autoregression, TV-VAR; 


量 超过 100 个 时 间 点 ), 且 目 前 只 能 用 于 两 变量 的 
建 模 分 析 (Piccirillo & Rodebaugh, 2019)。 感 兴趣 的 
研究 者 可 以 参阅 唐 文清 等 人 (2020) 以 及 Bringmann 
等 人 (2018) 的 文章 。 
4.5 密集 追踪 数据 分 析 方 法 的 选择 策略 

为 了 帮助 研究 者 针对 自身 研究 的 实际 情况 ， 
灵活 选择 适合 的 模型 方法 ， 这 里 对 各 种 主流 模型 
方法 的 优 缺 点 及 选择 策略 进行 归 类 和 总 结 ( 见 表 1 
及 图 6), 以 供 参 考 (Asparouhov et al., 2018; Gates 
et al., 2020; Gates & Molenaar, 2012; Piccirillo & 
Rodebaugh, 2019): 


5 讨论 


密集 追踪 方法 对 个 体 展开 多 个 密集 时 间 节 点 
的 数据 追踪 ， 能 更 好 地 研究 个 体 的 变化 过 程 和 内 
部 机 制 ， 在 外 部 效 度 、 数 据 分 析 方 法 等 方面 相 比 
传统 的 追踪 研究 具有 较 大 优势 。 为 了 鼓励 更 多 的 
研究 者 在 密集 追踪 数据 中 关注 动态 变化 的 过 程 、 


Bringmann et al., 2018)。TV-VAR 最 初 起 源 于 计量 
经 济 学 领域 , 但 近年 来 也 被 应 用 于 分 析 心 理 过 程 
的 动态 变化 (Bringmann et al., 2017; Bringmann 
et al., 2018; Chow et al., 2017; Haslbeck & Waldorp, 
2020)。 与 前 面 介绍 的 模型 方法 不 同 , TV-VAR 不 需 
要 假设 时 间 序 列 的 稳定 性 (Stationarity)， 而 是 允许 
模型 参数 随时 间 发 生变 化 ( 唐 文 清 等 ，2020; 
Bringmann et al., 2018)。 然 而 , TV-VAR 对 于 观测 
数据 点 的 数量 要 求 较 高 (一 般 需 要 对 每 个 被 试 测 


不 同时 间 节 点 间 存 在 的 自 回 归 效 应 等 ,本 文 主要 
介绍 了 两 种 密集 追踪 模型 方法 的 新 进展 ， 和 希望 能 
够 增进 研究 者 们 对 其 的 理解 、 运 用 和 进一步 发 展 。 
相 比 传统 的 分 析 方法 , DSEM 和 GIMME 都 更 多 地 
纳入 了 时 间 序 列 的 建 模 思想 ， 并 更 加 强调 对 于 个 
体 间 差异 的 关注 。DSEM 允许 研究 者 对 任意 数量 
的 随机 效应 进行 考虑 (Asparouhov et al., 2018)， 而 
GIMME 从 纯粹 的 基于 个 人 的 分 析 方 法 出 发 , 更 
直接 地 对 单个 个 体 的 模式 进行 探索 。 这 些 分 析 方 


表 1 密集 追踪 的 主流 方法 优 缺 点 


模型 方法 优点 


MLM/MSEM 简单 易 行 ; 对 时 间 点 数量 要 求 相对 较 低 


mlVAR 
动态 交互 ; 结果 可 视 性 强 

DSEM 可 以 将 随机 效应 分 解 为 个 体 和 时 间 两 个 方面 的 来 
源 ; 在 时 间 节 点 和 随机 效应 的 数量 上 限制 较 少 ; fo 
许 测量 时 间 节 点 之 间 间 隔 不 同 

(LV-)GIMME ”针对 单个 被 试 进 行 建 模 , 考虑 个 体内 部 的 动态 变化 
过 程 及 机 制 的 异 质 性 ; 适用 于 较 多 变量 间 的 动态 交 
互 ; 自动 化 搜索 


对 时 间 点 数量 要 求 相对 较 低 ; 适用 于 较 多 变量 间 的 假设 不 同 个 体 间 变量 的 相互 作用 机 制 是 同 


假设 不 同 个 体 间 变 量 的 相互 作 的 ; 难 
以 分 析 个 体内 部 的 动态 变化 过 程 及 机 制 ， 建 模 不 够 


时 
x 


考虑 潜 变量 因子 结构 


假设 不 同 个 体 间 变 量 的 相互 作用 机 制 是 
要 对 被 试 测量 较 多 时 间 点 


n 
z 
E 


可 能 出 现 模 型 过 拟 合 ; 无 法 将 几乎 不 随时 间 变 化 的 
变量 纳入 模型 ; 需要 对 单个 被 试 测量 较 多 时 间 点 


ik: GIMME: 组 迭代 多 模型 估计 (Group Iterative Multiple Model Estimation); LV-GIMME: 基于 潜 变量 进行 建 模 的 


GIMME (Latent Variable GIMME); MLM: 多 层 线性 模型 (Multilevel Model); MSEM: 多 层 结构 方程 模型 (Multilevel Structural 
Equational Model); mlVAR: 多 层 VAR 模型 (Multilevel VAR); DSEM : 动态 结构 方程 模型 (Dynamic Structural Equation Model). 
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图 6 密集 追踪 数据 方法 的 选择 策略 流程 图 


法 对 个 体 差 异 的 高 度 关 注 ， 反 映 了 近年 来 对 于 个 
体 特异 性 研究 的 倡导 和 发 展 的 趋势 。 


时 间 建 模 (Continuous Time Modeling) 可 以 较 好 地 
解决 时 间 离 散 化 假设 带 来 的 问题 , 但 它 需 要 涉及 


值得 注意 的 是 , GIMME 的 一 大 特点 在 于 其 模 
型 的 建立 是 探索 性 或 者 半 验 证 性 的 ,这 与 传统 心 
理学 研究 中 基于 前 人 研究 及 已 有 的 知识 经 验 提出 
假设 的 模型 、 利 用 大 规模 的 数据 收集 对 模型 进行 
进一步 验证 的 思想 存在 很 大 的 不 同 。 这 种 探索 性 
的 建 模 思 路 允许 研究 者 以 数据 驱动 的 方式 对 变量 
间 的 关系 进行 自由 探索 ,可 能 可 以 避免 先入为主 
的 固化 思维 , 为 研究 者 提供 更 多 新 的 发 现 和 启 
发 。 然 而 ,探索 性 分 析 常 常 也 伴随 模型 过 拟 合 、 
结果 难以 解释 等 问题 ,未 来 还 需要 研究 者 进一步 
探究 对 其 停止 搜索 规则 的 优化 设置 。DSEM 对 模 
型 中 随机 效应 的 数量 几乎 没有 限制 ,这 在 改善 模 
型 拟 合 的 同时 也 可 能 导致 模型 过 于 复杂 、 无 法 泛 
化 到 更 广泛 群体 的 问题 。 为 了 尽 可 能 减少 可 重复 
性 方面 的 问题 ， 建 议 研 究 者 在 实际 研究 前 对 其 实 
验 设 计 及 分 析 计 划 进 行 预 注册 ,并 在 文章 中 尽 可 
能 详尽 地 报告 其 数据 收集 、 处 理 及 分 析 等 过 程 。 

另外 ， 目 前 对 于 密集 追踪 数据 的 动态 分 析 大 
多 假设 测量 时 间 点 是 离散 化 的 ， 这 种 假设 在 连续 
时 间 测 量 ， 以 及 测量 时 间 间 隔 高 度 随 机 的 情况 下 
可 能 存在 一 定 的 问题 (Hamaker et al., 2018)。 连 续 


= 


复杂 的 微分 计算 (Voelkle et al.，2018)， 目 前 应 用 
研究 者 对 于 如 何在 实证 分 析 中 实现 连续 时 间 建 模 
还 不 够 熟悉 。 未 来 还 需要 更 多 研究 者 进一步 探索 
连续 时 间 模 型 的 简化 计算 与 自动 化 实现 ,为 应 用 
研究 者 提供 更 灵活 的 建 模 选 择 。 

密集 追踪 研究 关注 个 体内 行为 、 认 知 、 心 理 
等 的 动态 变化 过 程 , 这 些 内 在 过 程 及 机 制 往往 十 
分 复杂 ,可 能 存在 非 线 性 的 作用 模式 (如 ,两 变量 
间 存 在 交互 作用 )。 目 前 密集 追踪 数据 较为 主流 的 
分 析 方 法 大 多 仅 考虑 线性 效应 的 情况 (如 ，DSEM 
和 GIMME)， 但 仅 通过 线性 关系 对 这 些 动 态 效应 
进行 解释 很 可 能 是 不 够 的 。 最 近 , Kelava 和 Brandt 
(2019) 将 DSEM 与 动态 潜 类 别 分 析 (Dynamic Latent 
Class Analysis; Asparouhov et al., 2017) 相 结合 ， 提 
出 非 线 性 的 动态 潜 类 别 结构 方程 模型 (Nonlinear 
Dynamic Latent Class Structural Equation Modeling, 
NDLC-SEM)， 人 允许 在 组 间 和 组 内 水 平 进 行 非 线 性 
效应 的 分 析 (Kelava & Brandt, 2019)。 然 而 , 目前 
允许 对 密集 追踪 数据 进行 非 线性 建 模 的 模型 方法 
仍 较 为 有 限 ， 其 在 不 同情 境 下 的 具体 表现 还 不 明 
晰 ， 未 来 需要 更 多 方法 学 研究 者 对 其 进行 更 进 一 
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步 的 研究 。 

密集 追踪 研究 的 数据 收集 往往 基于 被 试 的 真 
实生 活 场 景 ,具有 较 高 的 外 部 效 度 。 然 而 , 真实 的 
测量 环境 往往 也 更 加 复杂 , 研究 者 可 能 需要 相应 
增加 对 潜在 干扰 因素 的 动态 测量 , 并 将 其 纳入 模 
型 进行 分 析 。 此 时 , 模型 中 涉及 的 变量 数目 较 多 
因此 采用 GIMME 对 数据 进行 分 析 可 能 是 更 合适 
的 。 然 而 , GIMME 对 于 单个 个 体 的 测量 时 间 点 数 
要 求 较 高 ， 这 大 大 增加 了 应 用 研究 者 密集 追踪 数 
据 的 收集 难度 。 未 来 可 以 进一步 探索 适用 于 高 维 
密集 追踪 数据 、 且 对 测量 时 间 点 要 求 相 对 较 低 的 
模型 方法 ， 以 及 进一步 发 展 密集 追踪 数据 的 收集 
方法 ， 尽 可 能 降低 其 收集 的 难度 。 

总 之 ,密集 追踪 方法 的 一 大 优势 在 于 探究 个 
体内 部 的 动态 变化 过 程 及 其 机 制 。 虽 然 目 前 密集 
追踪 数据 分 析 在 方法 学 领域 已 有 较 多 新 的 进展 ， 
但 它 在 实际 研究 中 应 用 不 足 , 许多 应 用 研究 者 缺 
乏 对 于 时 间 这 一 要 素 的 关注 (Cronin & Vancouver, 
2019; Voelkle et al., 2018; Zhou et al., 2020)。 现 在， 
较 多 密集 追踪 数据 的 分 析 方 法 已 经 可 以 方便 地 通 
过 软件 实现 (如 , DSEM 和 GIMME), 大 大 减少 了 
实证 研究 应 用 的 阻碍 。 和 希望 本 文 可 以 帮助 实证 研 
究 者 们 了 解 、 掌 握 更 灵活 的 建 模 和 分 析 方 法 ,更 
好 地 发 挥 密集 追踪 数据 的 优势 ， 尽 可 能 减少 对 数 


据 结果 的 误 读 。 
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Intensive longitudinal data analysis: Models and application 


ZHENG Shufang, ZHANG Lijin, QIAO Xinyu, PAN Junhao 
(Department of Psychology, Sun Yat-sen University, Guangzhou 510006, China) 


Abstract: In the fields of psychology, education, and clinical science, researchers have devoted increased 


attention to the dynamic changes and personalized modeling of individuals' behaviors, minds, and treatment 


effects over time. Intensive longitudinal data is a set of measures collected at multiple time points with 


higher frequency over shorter periods. Thus, it can be used in the analysis of the dynamics and mechanisms 


of within-person processes. In recent years, intensive longitudinal design has become one of the most 


prominent and promising approaches in psychological research. However, many of these researches still rely 


on traditional data analysis methods. Many models have been proposed to analyze intensive longitudinal 


data, including top-down approaches (e.g., dynamic structural equation model, DSEM) and bottom-up 


approaches (e.g., group iterative multiple model estimation, GIMME). Both of the methods can 


conveniently model autoregressive and cross-lagged effects in intensive longitudinal data. 


Key words: intensive longitudinal data, time-series, DSEM, GIMME 
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附录 A: MLM 5 RDSEM 
Mplus 代码 


1. MLM 代码 : 
DATA: 
FILE = [DATA LOCATION PATH] two-level data 
with Trend.csv; 
VARIABLE: 
NAMES ARE URGE DEP JS HS PERSON TIME; 
USEVARIABLES URGE DEP JS HS TIME; 
cluster-PERSON; 
BETWEEN - JS HS; 
WITHIN = TIME DEP;! Time is only used in 
the Within-Level model; 


DEFINE: 
CENTER DEP (groupmean); 


ANALYSIS: 
TYPE = TWOLEVEL RANDOM;!Two-Level 
Model with paths that are latent variables; 


MODEL: 
%WITHIN% 
bli | URGE ON DEP; 
! urge to smoke is regressed on Lag-l urge to 
smoke 
trend | URGE ON TIME;! urge is regressed on 


time to capture the trend over time, the !slope is latent; 


%BETWEEN% 
URGE bli; 
URGE with bli TREND;bli with TREND; 


[URGE]; !mean intercept of urge to smoke; 
[bli]; !mean of the regressive slope; 


[TREND]; 
URGE; !urge to smoke intercept variance; 
bli; !regressive slope variance; 


TREND; 


URGE on JS HS; 


!The person-specific intercept is predicted by 
Job Stress and Home Stress; 

bli on JS HS; 

!The regressive slope is predicted by Job Stress 
and Home Stress; 


TREND on JS HS; 


OUTPUT: TECH! TECH8; 


2. RDSEM 代码 : 
DATA: 
FILE = [DATA LOCATION PATH] two-level data 
with Trend.csv; 
VARIABLE: 
NAMES ARE URGE DEP JS HS PERSON 
TIME; 
USEVARIABLES URGE DEP JS HS TIME; 
Lagged = URGE (1) DEP (1);!Created a Lag-1 
variable for urge to smoke and depression; 
CLUSTER=PERSON; 
BETWEEN = JS HS;! Job Stress and Home 
Stress are only used in the Between-Level model; 
WITHIN = TIME;! Time is only used in the 
Within-Level model; 


ANALYSIS: 

TYPE = TWOLEVEL RANDOM; !Two-Level 
Model with paths that are latent variables; 

ESTIMATOR = BAYES;!Lagged variables can 
only be estimated with Bayes in Mplus; 

BITERATIONS = (1000);!Run at least 1000 
iteration of the MCMC algorithm; 

BSEED=1028; 


MODEL: 

%WITHIN% 

PHI | URGE^ ON URGE“ 1; 

! urge to smoke is regressed on Lag-l urge to 
smoke 

PHI2 | DEP“ ON DEP“1; 

! depression is regressed on Lag-1 depression 

PHI3 | DEP^ ON URGE“; 


! Depression is regressed on Lag-1 urge to 
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PHI4 | URGE^ ON DEP^I; 

! urge to smoke is regress on Lag-1 depression 

TRENDI | URGE ON TIME;! urge is regressed 
on time to capture the trend over time 

TREND2 | DEP ON TIME:! dep is regressed on 
time to capture the trend over time 

LogV1 | Urge; 

! the Within-Level residual variance of urge to 
smoke 

LogV2 | DEP; 

! the  Within-Level 


residual variance of 


depression 


%BETWEEN% 

URGE]; !mean intercept of urge to smoke 
DEP]; !mean intercept of depression 

[PHI1]; !mean of urge autoregressive slope 
PHI2]; !mean of depression autoregressive 
slope 
PHI3]; !mean of dep on urge&1 slope 
[PHI4]; !mean of urge on dep&1 slope 
TREND1]; 

TREND2]; 

LogV1]; 

‘log of the mean of the urge residual variance 


LogV2]; 


!log of the mean of the urge residual variance 


URGE; !urge to smoke intercept variance 

DEP; !depression intercept variance 

PHI1; !urge to smoke autoregressive slope 
variance 

PHI2; !depression autoregressive slope variance 

PHI3; !dep on urge&1 slope variance 

PHI4; !urge on dep&l slope 

TRENDI; 

TREND2; 

LogV1; 

!between-person variance of urge to smoke 
residual variance 

LogV2; !between-person variance of depression 


residual variance 


URGE on JS HS; 

DEP on JS HS; 

!The person-specific intercepts are predicted by 
Job Stress and Home Stress; 

PHII on JS HS; 

PHI2 on JS HS; 

PHI3 on JS HS; 

PHI4 on JS HS; 

!The person-specific slopes are predicted by Job 
Stress and Home Stress; 

LOGVI on JS HS; 

LOGV2 on JS HS; 

!The person-specific residual variances are 


predicted by Job Stress and Home Stress; 


URGE DEP PHI1-PHI4 LOGVI LOGV2 WITH 
URGE DEP PHI1-PHI4 LOGV1 LOGV2; 


OUTPUT: TECHI TECH8 STDYX; 


附录 B: GIMME 实证 演示 R 代码 


library (DiagrammeR) 
library (gimme) 
library (mvtnorm) 
library (dplyr) 

library (xlsx) 


#1、 数 据 预 处 理 

wd<-"~"# 设置 文件 夹 路 径 
setwd (wd) 

dat<- read.csv ("Two-Level 
header-F) 4 读 取 数据 文件 
colnames (dat)<-c ("URGE", "DEP", "JS", "HS", 
"PERSON", "TIME") # 指定 各 变量 名 ,其 中 URGE 指 
吸烟 欲望 , DEP 指 抑郁 水 平 

datawd<-paste0(wd, "/data") 


Data with  Trend.csv", 


dir.create(datawd) 

setwd (datawd) 

for (i in 1:length (unique (dat$PERSON)))1 

8 选取 个 体 i 的 数据 进行 预 处 理 

person dat--dat[which (dat$SPERSON--i), c ("URGE", 
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"DEP", "TIME")] 


# 分 别 去 除 变 量 URGE 和 DEP 随时 间 变 化 的 整体 趋势 


(detrend) 

temp = lm (person dat$URGE~person dat$TIME) 
person_dat$urge<-residuals (temp) 

temp2 = Im (person dat$DEP-person dat$TIME) 


person dat$dep--residuals (temp2) 


# 将 个 体 i 的 数据 单独 写 入 excel 文档 "data i.csv" 
write.csv (person dat[, c ("urge", "dep")], 


paste ("data", i, ".csv"), row.names = FALSE) 


#2、 利 用 gimmeSEM pk Bc h PEAS IS HEC £s TE 
fit<-gimmeSEM ( 
data = datawd, # 存放 数据 的 文件 夹 路 径 


out = paste0(wd, "/output"), # 存放 输出 结果 的 文件 
夹 路 径 

sep=",", # 声明 数据 分 隔 符 , 这 里 csv 格式 的 数据 
分 隔 符 为 逗号 

header - T, 4 声明 数据 中 是 否 包括 变量 名 

ar = TRUE, # 指定 模型 中 是 否 包 括 自 回归 ， 默 认为 
TRUE 

plot = TRUE, s 指定 结果 是 否 输出 模型 图 ， 默认 为 
TRUE 

groupcutoff —75 4 路 径 可 以 被 纳入 群体 模型 的 临界 
值 (显著 改善 个 体 模型 的 最 小 比例 ， 默 认为 75%) 


) 
fit2<-aggSEM( 


data = datawd, # 存放 数据 的 文件 夹 路 径 
out = paste0(wd,'/output'), # 存放 输出 结果 的 文件 夹 路 径 
sep = 


header = T) 


